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slope  of  curve  of  pitching  aoaent  coefficient  vs  oc 
change  in  per  degree  pitch  control  incidence,  6p 

noraal  force  coefficient 

slope  of  curve  noraal  force  coefficient  C^vs  o< 
change  in  per  degree  pitch  control  incident, 
yawing  aoaent  coefficient 

slope  of  curve  of  yawing  aoaent  coefficient,  vs  p 
change  in  per  degree  yaw  control  incidence,  6y 
change  in  C,^  per  degree  roll  control  incidence, 
side  force  coefficient 

slope  of  curve  of  side  force  coefficient  vs  ^ 

change  in  per  degree  yaw  control  incidence , 

change  in  Cy  per  degree  roll  control  incidence, 

reference  length  for  coefficients  «  2ft 

aoaent  of  inertia  about  T-  axis 

B 
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kyp 


f 

ft 

r 

*e 

p 

S 


H 

V 


V 

V 
w 

Y« 

Yv 


■oaent  of  inertia  about  axis 

■oaent  of  inertia  about  7^  axis 

CBTT  autopilot  coordinator  branch  gain 

roll  rate  about  7^ 

roll  acceleration  about  7^ 

constant  or  eguilibriue  roll  angular  rate 

dynaelc  pressure 

pitch  rate  about  j 

pitch  angular  acceleration  about  y 

constant  or  eqnilibrina  pitch  angular  rate 

yaw  angular  rate  about  z^ 

yaw  angular  rate  coesand  (coordination  cossand) 
yaw  angular  acceleration  about  7g 
reference  area  for  coefficients  *T  ft* 

▼eloclty  coeponent  in  Indirection 

▼elocity  coeponent  in  y  direction r  assueed  to  be  constant 
constant  eissile  flight  path  relocity 
aissile  welocity  vector 
velocity  coeponent  in  z  direction 

body-fixed  roll  axis,  along  axis  of  syeeetry,  positive  forward 

body-fixed  pitch  axis,  positive  starboard 

vehicle  axis  in  local  horizontal  direction,  approxiaated  as 
inertial  axis 

body- fixed  yaw  axis 
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▼•hid*  axis  in  detmvard  directiaa  along  local  gravity 
vactor*  approxiaatad  as  inertial  axis 


'nz 

\ 

4c 

4 

4e 

e 

4» 

^Pc 

^Yc 


achievad  noraal  acceleration  in  z_ direction 
coaaanded  noraal  acceleration  in  ^direction 

O 

achievad  noraal  acceleration  in  7  direction 

s 

achievad  noraal  acceleration  in  z*y direction 

achieved  noraal  acceleration  in  y  direction 

noraal  acceleration  coaaand  froa  gaidance  coaputer 
in  Zydlrectlon  plas  anti -gravity  bias  coaaand 

noraal  accelention  gaidance  coaeand  in  Zydirection 

noraal  acceleration  gaidance  coaaand  in  indirection 

roll  attltade  coaaand  froa  gaidance  coapater#  zero  degrees 
in  -z  direction  and  90  degrees  in  indirection 

roll  attltade,  zero  degrees  in  -z^ direction  and  90  degrees 
in  y^  direction 

roll  attltade  error. 

Elevation  Enlar  tnglar,  aacond  rotation p  ~  rsln  ^  )dt 

iziaath  Baler  ingle,  first  rotation 
aboat  Yy  f  J {^ain  ^  ♦  rcos  p  )  dt 

pitch  control  incidence  (positive  tail  incidence  prodaces 
negative  pitiching  aoaent) 

coaaanded  pitch  control  Incidence, 

yaw  control  incidence  (positive  tail  incidence  prodaces 
negative  yawing  noaenta) 

coaaanded  yaw  control  incidence, 

roll  control  Incidence  (positive  tail  incidence  prodaces 
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posit iT«  rolling  lossnt) 
coaaanded  roll  control  Incidenca 


constant  or  eqailibrlua  angle-of* attack 

angle-of-attack 

angle-of-attack  rate 

Bodlfled  fora  of  estimated  angle*of-attack  for  autopilot 
coordination  coaaand 

angle  of  sideslip 

sideslip  angular  rate 


! 

I 
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I.  EABiaagBs  ssasiimii  iJiimis  as  4  siisriszEm  aiasiLi 

The  deterainatlon  of  changes  in  systee  perforaance  due 
to  changes  In  pacaneters  is  of  great  Isportance  in  engi¬ 
neering  analysis  and  design.  Sensitivity  questions  arise 
when  eodel  uncertainty  is  present,  or  when  a  range  of  oper¬ 
ating  conditions  is  contesplated. 

The  questions  of  paraeeters  sensitivity  particularly  arise 
in  the  fields  of  engineering  whera  aodels  are  idealized, 
inexactly  identified,  or  the  systeas  theaselves  are  subject 
to  unpredictable  changes  with  tiae  due  to  envlroaental, 
aaterial  property  or  operational  influences  so  that  there  is 
always  a  discrepancy  between  the  physical  reality  and  the 
aatheaatlcal  aodel. 

Sensitivity  analysis  provides  the  angineer  with  aethods  for 
investigating  or  ainlaizing  such  paraaeter  deviations. 

In  general,  the  diagrae  of  a  systea  can  be  represented 
by  a  single  block  as  given  on  Pig.  1.1. 


SYSTEM 

Figure  1.1  General  Representation  of  a  Dynaalc  Systea. 


Proa  a  aatheaatlcal  point  of  view,  what  we  call  a 
systea  Is  the  expliclty  or  iapliclty  given  relationship 
between  the  input  signal  u  and  output  signal  y.  In  general, 
u  and  y  can  be  vectors.  The  character  of  this  relationship 
is  coanonly  called  the  structure  of  the  systea. 
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For  Qxaapler  the  stractare  of  the  systen  may  be 
characterized  by  the  order  of  a  differential  or  difference 
equation,  linearity  or  nonlinearity,  the  order  of  the  numer¬ 
ator  and  denominator  of  a  rational  transfer  function  and  the 
quantitative  properties  of  the  system  parameters. 

Typical  parameters  are  initial  conditions,  time- 
invariant  or  time-variant  coefficients,  natural  frequencies 
and  sampling  periods. 

The  change  of  the  state  or  the  change  of  the  output 
variable  with  time,  can  be  caused  by:  (1)  the  influence  of 
input  signals,  (2)  the  change  of  parameters.  These  quanti¬ 
ties  are  shown  in  Fig.  1.2. 


Figure  1.2  Quantities  Affecting  the  Dynamics  of  a  System. 

This  work  addresses  the  application  of  the  sensitivity 
analysis  to  the  variation  of  aerodynamic  paraaeters  of  a 
Bank-to-Turn  missile. 

PARABETEB  SENSITiyiTT  is  the  effect  of  parameter  changes  on 
the  dynamics  of  a  system,  say,  the  time  response,  the  state, 
or  any  other  quantity  characterizing  the  system  dynamics. 

In  order  to  have  a  realistic  model  to  work  with  a  RASA 
Contractor  Report  (Ref.  2]  was  adopted  as  a  reference  for 
application  of  the  sensitivity  analysis  on  the  variation  of 
aerodynamic  parameters  of  a  Bank-to-Turn  missile.  Appendix 
A,  B  and  C  give  the  detailed  aerodynamic  assumption  of  the 
system. 


One  brief  explanation  of  the  sensitixity  theory  and 
analysis  is  given  as  a  guideline  to  better  unierstand  the 
snbsegaent  work. 

&n  analysis  of  the  sensitivity  of  the  aerodynanic  paraae- 
ters  applied  to  the  case  of  the  linear  uncoupled  pitch  and 
coupled  roll-yaw  autopilot  is  performed. 

The  pitch  and  roll-yaw  autopilots  are  coupled  as  given 
in  Flg.C.1  in  the  appendix  C.  &  nonlinear  parameter  sensi¬ 
tivity  analysis  is  performed. 

Conclusions  about  the  influence  of  parameters  of  concern 
are  given  and  a  comparison  between  the  linear  and  nonlinear 
case  is  done.  Comments  and  recommendations  for  the  future 
are  outlined  in  order  to  delineate  the  continuity  of  the 
present  work. 
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iz.  sigsmmi  iijS2£i 


i.  1BTBOO0CTIOH 

The  basis  of  all  sensitivity  considerations  in  the  case 
of  tise-invariant  paraseter  variations  is  ^ha  so-called 
sensitivity  fnnctions.  Dynaaic  systess  can  be  characterized 
in  several  ways: in  the  tiae  doaain,  in  the  fregaency  doaain, 
or  in  terns  of  a  perfoaance  index.  There  is  evidently  an 
adagnate  naaber  of  ways  to  define  the  sensitivity  functions 
of  a  dynaaic  systea.  The  definition  that  is  actually  used 
depends  on  the  fora  of  the  aatheaatical  aodel  as  well  as  on 
the  purpose  under  consideration. 

The  sensitivity  functions  can  be  classified  into  the 
following  three  categories: 

(1)  Sensitivity  functions  in  the  tiae  doaain 

(2)  Sensitivity  functions  in  the  freguency  doaain,  and 

(3)  perfoaance- index  sensitivity 

In  this  chapter  we  will  outline  just  the  sensitivity  func¬ 
tions  in  the  tiae-doaain  of  the  continuous  systeas. 
Iplications  of  this  analysis  follow  throughout  this  work. 

The  reader  can  find  detailed  inforaations  about  the 
others  two  categories  of  sensitivity  functions  in  the 
CBef.  1]. 

B.  TBiJBCTOBT  SBISITITITT  POBCTIOB  OF  COBTIBOOOS  STSTBBS 

k  continuous,  possibly  nonlinear  systea  of  nth  order 
can,  in  general,  be  described  in  the  state  space  by  a  vector 
differential  equation  as  seen  in  Ego. 2.1. 
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X  -  ,  X(t^)  »  XO  (2.1) 

Here  X  is  an  n  x  1  state  xector,  f  an  n  x  1  sector  func¬ 
tion*  ?  an  input  vector*  a  noainal  r  x  1  parameter 
vector*  and  X  is  the  n  x  1  initial  condition  vector  or 
initial  state.  Bgn.2.1  is  called  the  NOHIN&L  STATE 
EQUATION . 

Assueing  that  the  parameter  vector  deviates  from  the 
nominal  value  by  ^o^*  we  have  the  Egn.2.2. 


X 


f(X*t,a*«c),  x(to) 


XO 


(2.2) 


with  the  initial  conditions  XO  anchanged.  This  equation 
is  called  the  ACTUAL  STATE  EQUATION. 

Now  it  is  assumed  that  Eqn.2.2  has  a  unique  solution,  X» 
Z(t*^)  for  all  admissible  initial  conditions  and  parameter 
values. 

X  is  of  coarse  a  function  of  u  ,  x®  and  t  as  well. 

O 

However  *  this  dependence  is  not  needed  for  the  following 
considerations  and  will*  therefore,  be  dropped  for  ease  of 
notation.  Furthermore*  the  solution  X  is  assumed  to  be  a 
bounded  continuous  function  in  t  a:.d  ^  . 

If  the  parameter  takes  on  its  NOHINAL  value  »  the 

noainal  solution  X^  »  X(t*oJ^)  is  obtained.  If,  on  the 
other  hand  the  ACTUAL  solution  is  given  by  X  ^  X(t,  *<  )  , 
then  the  para me ter- induced  change  of  the  state  vector  is 
given  by 


A  X(t*o<)  «  X(t*  •<)  -  X(t*  o^) 


(2.3) 
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A  f irst-ocder  approxiaation  of  a  X  can  be  obtained  by 
using  Taylor  expansion  in  the  form  of  Eqn.2.4. 


AX(t,‘5c)  * 


y 

j = 1 


(2.4) 


This  equation  can  be  viewed  as  a  definition  of  the 
parameter-induced  trajectory  deviation. 

The  subscript  oc©  shall  indicate  that  the  partial  deriv¬ 
ative  expressed  by  *7)  is  taken  at  nominal  parameter  values. 

Let  the  state  X  of  a  continuous  system  be  a  continuous 
function  of  a  time- invariant  parameter  vactor  o<,  = 
{ oc,  ,  otj  .  Then  the  partial  derivative  can  be 

defined  as 


-  — 


(2.5) 


i®) f  2f  .  .  .  #  n 
j*l,  2,  .  .  .  ,  r 

#qn.2.5is  called  the  trajectory  sensitivity  vector  with 
respect  to  the  jth  parameter. 

Hote  that  the  trajectory  sensitivity  vector  is  of  the 
same  dimension  as  the  state  vectorr  namely,  n.  Its  compo¬ 
nents  are  the  TRAJECTORY  SEHSITIYITY  FONCTIONS  as 


7  -  A'DXi(t,5} 


(2.6) 


Bqn.2. 6  is  the  partial  derivative  of  the  ith  state  vari* 
able  in  relation  to  the  jth  parametar  as 


I 


Vo 


(2.7) 


Hence,  all  n  x  r  trajectory  sansitixity  functions  fora 
the  trajectory  sensitivity  aatrix  as  given  in  Eqn.2. 8  or 


^  ^  ^  X 

A  *  {  A,  ,  A  )  =  -ir--—  I  — 

^  Q  oc  <»<>- 


^  Xj 
D«»<v 


(2.8) 


(2.9) 


^Xn  ^Xn 


The  columns  of  A  are  the  trajectory  sensitivity 
vectors  .  Here  A  is  the  Jacobian  aatrix  of  the  state 
vector  with  respect  to  the  paraaetar  vector  <K  ,  taken  at 
the  noainal  paraaeter  values. 

With  these  definitions  the  paraaeter-indaced  change  of 
the  trajectory  can  be  taken  as 


X(t,oC)  «  A  (t,os,) 


(2.  10) 


Where  o<  ■  ^  ,  which  is  the  kCTOiL  paraaeter 

vector  of  the  systea. 


C.  TBAJECTOBT  SEBSITIVITT  EQOITZOVS  OP  COMTIBOOJS  STSTEBS 

Lets  consider  the  general  contlnaoas  system  described  as 
previously  by  the  state  equation  (Bqn.2.1). 


X  »  f(X,t,u,  OC),  X(t^)  »  XO 


(2.11) 


Where  X  denotes  the  n-dioensional  state  vector/  o<  the 
r-diaensional  parameter  vector,  £  an  n-dimensional  vector 
function,  and  u  the  input  vector  independent  of  oC  .  It  is 
assumed  that  the  continuity  .conditions  are  fulfilled  and 
that  ot,  is  time-invariant. 

Considering  o<  -parameters  and  taking  the  partial 
derivative  of  X  (Eqn.2.11)  with  respect  to  ,  one  obtains, 
by  the  application  of  the  chain  rule 

- -  —  - - ♦ - ^  - - a  0  (2. 12) 

"b  1)  X  ^  J  ^  cx^' 

The  derivative  of  the  initial  conditions  vector  X®  with 
respect  to  is  zero,  since  X®  does  not  depend  on  o<  . 

If  oc  is  r-diaensional,  there  are  r  equations  of  the 
form  of  Eqn.2.12.  If  we  now  interchange  the  sequence  of 
taking  the  derivative  with  respect  to  tine  t  and  ,  and 
then  let  o<.  approach  to  ,  one  obtains  the  following 

equation 


X 


:u* 

'j 


0 


(2.13) 
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Here  Aj  trajectory  sensitivity  vector 

with  respect  to  the  jth  paraeeter. 

The  Eqn.2. 13  is  called  the  sensitivity  eqaation  in  the 
state  space  or  the  TRAJECTOHI  SESSITITITY  EQOITION. 

The  above  shows  that  for  o<.  -paraaeters  all  initial 
conditions  of  the  trajectory  sensitivity  eqnations  are  equal 
to  zero. 

How  consider  the  output  vector  equation,  as  given  by 


g(X,t,a,«) 


(2.14) 


In  a  procedure  siailar  to  that  above  one  obtains  the 
algebraic  sensitivity  equation  as  seen  in  Eqn.2. 15. 


\J  I  a  I  ^  , 


♦  • 


(2.15) 


output  sensitivity  vector  » 


i'  -  This 


^  _  Which  relates  the 
ON 

to  the  trajectory  sensitivity  vector 
equation  is  called  the  YECTOB  OHTPar  SEHSITIYITY  EQUATION. 

Using  the  trajectory  sensitivity  matrix  A  and  the 
output  sensitivity  aatriz  (r*  ,  the  above  result  can  be 
rewritten  in  the  following  general  fora 


A  »  — I-  ^  1_ 


(2.16) 


<r  .  3L.,  X 
7)5  v. 


2L,  _ 

*7)  o<. 


(2.17) 
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These  eqaations  are  called  the  STATE  SESSITIVITY 
BQOATIONS  of  the  systee.  It  is  seen  that  these  equations 
are  linear  whether  the  original  systee  is  linear  or  nonlin¬ 
ear. 

If,  in  particular,  the  original  systee  is  linear,  the 
state  equations,  take  the  fore 


X  a  A  X  ♦  B  u  ,  X(to)  *  XO 


(2.18) 


r  «  c  X  ♦  D  a 


(2.19) 


ihere,  in  general, 

A»A(oC)  ,  b'  *  B(«.)  ,C«C(o<)  ,D»d'(<), 

X  =  X(t,CK  )  ,  aid  7  a  7(t,ix  )  . 

Note,  however,  that  7  is  not  a  function  of  oc  if  u  is 
defined  as  an  external  input  of  the  systea. 

Now  taking  the  partial  derivatives  with  respact  to  ©Sj  » 
reversing  the  order  of  differentiations  with  respect  to 
tiae  and  ofj  ,  and  letting  o?  approach  to  ,  the 
TSAJECTORT  SBNSITIVITT  EQOATIONS  are  obtained 


7  X  .  7)  A 


-  Db  - 

LI®  ♦  Oft) 

«Ro  ®<o 


.  ACto)  =  0  (2.20) 


Where  A^  a  KoC/,}  ,  X^  ■  X(t,  o^”^)  ,  and  j»  1,  2,...,  r. 
The  initial  condition  vector  XjV.‘C:e)»  is  again  zero  since 
T(to)  does  not  depend  on  ^  • 
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By  slallar  pcocadora  the  VBCT3B  SBNSITI7Iir  BQO&TIOH 
becoees 


1.L 


LI* 

«^o 


(2.21) 


Where  «  C(o<o)  c. 

In  the  case  of  linear  systees,  the  vector  sensitivity 
equations  have  the  saee  i  eatriz  as  the  noeinal  state  equa¬ 
tions  and  hence  the  saee  characteristic  eatriz  si  -  &. 

They  differ  froa  the  noeinal  original  state  equation  only 
in  the  driving  function  and  the  Initial  conditions.  The 
latter  are  all  zero.  The  driving  functions  can  be  obtained 
by  solving  the  noeinal  state  equations. 

Here,  in  the  linear  case,  the  sensitivity  equations  can 
be  joined  to  the  systee  equations  ,  forelng  the  so-called 
CX3HBI1IBD  STSTBH  as  given  in  Bqns.2.22  and  2.23. 


(2.22) 


(2.23) 


The  sleultaneoue  solution  of  these  differential  equa¬ 
tions  using  the  saee  standard  procedure  for  each  of  the 


above  aatriz  egaatlona  yields  the  oatpat  vector,  the  state 
vector,  and  the  corresponding  oatpat  and  trajectory  sensi¬ 
tivity  vectors  ,  and  Aj  •  If  there  are  r  paraaeter 
variations,  r  systeas  of  agnations  of  the  above  type  have  to 
be  solved.  Since  the  hoaogeneous  part  of  the  original 
differential  agnation  is  identical  with  the  hoaogeneons  part 
of  the  sensitivity  egoations  with  respect  to  all  paraaeters. 

1  graphical  interpretation  of  Egns.2.22  tnd  2.23  is 
given  in  Plg.2.1. 

Since  the  driving  function  of  the  sensitivity  aodel 
contains  the  noainal  state,  the  aeasnring  circuit  for  the 
trajectory  sensitivity  functions  consists  of  a  connection  of 
the  noainal  original  systea  and  the  sensitivity  aodel  as 
Illustrated  by  Fig. 2.1. 

If  the  actual  input  *u  is  applied  to  this  stricture  the 
trajectory  and  output  sensitivity  vectors  Xj  and  (Tj  can  be 
aeasured  at  the  points  1  and  2  respectively  in  the  Fig. 2.1. 

In  order  to  aeasure  all  r  sensitivity  vectors  siaultane- 
ously,  r  sensitivity  aodels  are  reguirad. 

D.  STBOCTOBiL  HBIHOO 

The  basis  for  the  aeasureaent  of  the  trajectory  sensi¬ 
tivity  functions  is  the  STBDCTDBAL  IHTEBPBETillON  of  the 
trajectory  sensitivity  eguation. 

The  physical  systea  represented  by  the  trajectory 
sensitivity  egoations  is  called  the  trajectory  sensitivity 
aodel  or  the  sensitivity  aodel  in  the  stare  space. 

Begardless  of  the  nature  of  the  original  systea  the 
sensitivity  aodel  in  the  state  space  is  always  linear.  The 
graphical  illustration  of  sensitivity  eguation  can  be  given 
by  Fig. 2. 2.  The  systea  described  by  the  sensitivity  eguation 
is  refered  to  as  the  sensitivity  aodel  of  the  original 
systea.  For  each  output,  a  systea  has  as  aany  sensitivity 
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■odals  as  parassters  of  Interest.  Both  the  nosinal  original 
systes  and  the  corresponding  sensitivity  sodels  fors  the 
OOHBIHED  STSTBB.  The  sensitivity  solel  is  always  linear.  If 
T(t««(w)  is  the  ontput  of  a  systep^  the  corresponding  sens!* 
tivity  functions  (tyoC^)  ■  outputs  of 
the  corresponding  sensitivity  sodels.  Consequently,  in  order 
to  seasure  the  output  sensitivity  functions  sisultane- 
ously,the  nosinal  original  systes  and  the  sensitivity  sodels 
have  to  be  seasured  at  the  outputs  of  the  sensitivity 
sodels. 

The  double  frase  used  for  the  original  systei  in  Fig.  2.2 
is  to  indicate  that,  in  general,  the  systes  say  be  nonlin¬ 
ear,  whereas  the  sensitivity  sodels  are  always  linear.  In 
the  nonlinear  systess  the  Eqn. 2.12  applies  and  different 
prograss  have  to  be  set  up  for  the  original  and  the  sensi¬ 
tivity  equations. 

For  sose  applications,  the  sensitivity  functions  with 
respect  to  all  paraseters  are  required  slnultaneously.  If 
there  are  r  paraseters  of  Interest,  r  sensitivity  sodels 
would  be  needed.  This  isplies  a  rather  extensive  cosputer 
tise. 

A  sethod  of  deterslnlng  all  output  sensitivity  functions 
of  a  systes  sisultaneously  by  a  single  sensitivity  sodel  is 
available  ,  which  is  also  called  the  sethod  of  SENSITIFITT 
POINTS.  This  sethod  has  application  just  for  linear  systess. 
Detailed  explanation  about  the  application  of  the  sensi¬ 
tivity  points  theory  can  be  found  in  the  [Bef.  1]. 

B.  OOTPUT  SEISITITITT  FONCTIOI  OF  CONTIIOOOS  STS TEES 

Consider  the  input-output  behavior  of  a  continuous, 
possibly  nonlinear,  single-variable  systes  described  by  an 
ordinary  differential  equation  of  the  type 
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f{  »y*"  0<^  )  ■  0  (2.24) 

with  th«  initial  condition  y{t^)  »  y.  ®»  i»0,  1,...n-l.  y 
danotas  tha  ootpnt  signal#  t  the  tiae#  and  o<  a  single 
tisa* invariant  or  slowly  varying  -paraaeter  that  has 
noainal  value  «<«  . 

In  general#  f  is  a  function  of  the  input  u  as  well. 
Bowavar#  if  u  is  an  external  input  which  does  not  depend 
on  #  tha  dapandenca  of  f  on  n  is  not  relevant  in  further 
considerations. 

Let  ns  suppose  that  tha  above  differential  equa¬ 

tion  has  the  unique  solution  given  below 


To  «  T(t#  o<o)  (2.25) 

Which  one  shall  call  the  NOMINAL  SOLUTION. 

Let  us  now  assuae  that  the  paraaater  changes 
froa  to  X*  ♦  ^®^#  where  is  tiae-invariant  or 

slowly  varying  with  tine.  o<,  is  called  the  actual  paraneter 
value.  The  corresponding  ACTUAL  DIFFERENTIAL  EQUATION  can# 
then#  be  written  as 

(n)  (n-i; 

f  ■  {  J)  #y  ^•••«y«t#o(}  *  0  (2.26) 


Note  that  by  this  change  of  into  oc  the  Initial 
conditions  renain  unchanged#  naaely  y(t^)  *  y.<>. 

^  L.* 
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Th«  correspondi&g  solution  is  given  as 


y  «  y(t,o<) 


(2.27) 


Which  ve  shall  call  the  iCrOAL  (or  PBBT05B2D}  SOLOTIOH. 

It  is  assased  that  y(t,  «<)  is  of  the  saae  type  as 
y(t«^)»  and  y(t,  o<)  deviates  infinitesiaally  fcos  y(t,o(o) 
if  ^  deviates  infinitesiaally  froa  •  The  conditions 
for  fulfilling  this  reguireaent  are  given  in  the  aatheaat* 
ical  literature.  For  our  purpose,  it  is  sufficient  to  know 
that  y  is  continuous  in  ^  if  f  is  continuous  in  y  which  is 
true  for  all  continuous  systeas  and  t  0. 

With  the  above  assnaptions,  tha  actual  solution  7  (t,  ocg 
*Zot)  can  be  expanded  into  a  Taylor  series  around  , 
yielding  the  Bgn.  2.28. 


y  (t  ,«c) 


y(t,^  ♦  -i-|Aoc4  -1- 


(2.28) 


If  lkot«  oCq  ,  the  Taylor  series  can  be  truncated  at  the 
linear  tera. 

This  gives  the  Bgn. 2.29. 


y  (t/w) 


y  (t  ,o(o) 


♦ 


I  ^  «< 


(2.29) 


For  finite  values  of  A  n<  ,  this  expression  can  be 
considered  a  first-order  approxiaation  of  y(t,o(). 

The  actual  output  can  be  written  as 
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y(t,oi)  t  7(t,(Ko)  ♦  O'  (t,o<o)  AOC 


(2.30) 


Where  y(t»^e)  la  the  noelnal  oatpat  and  (T  (t,  <i<^)  is 
the  paraeeter* induced  output  error. 

With  0“  (tf  Xo)  defined  paraeeter 

Induced  output  error  is  in  these  teres  defined  as 

y(t,<»<)  t  O’ct.oto)  Ao<.  (2.31) 
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III.  0£  Ssssisillll  MILISIS  £2  UMSiS  stitehs 

i.  IITBODOCTION 

The  sensitivity  theory  described  previously  will  be 
applied  to  the  case  of  an  uncoupled  pitch  autopilot  and  to 
a  roll-yaw  coupled  autopilot  presented  in  the  appendix  C. 
For  each  case,  the  linear  equations  of  the  noeinal  system 
are  presented  in  state  variable  fore  and  the  correspondent 
trajectory  sensitivity  equations  are,  derived.  A  sensitivity 
analysis  of  the  systees  is  performed. 

B.  OICOOPLED  PITCH  AOTOPILOT  AMALTSIS 

1,  Linear  Equations  of  the  Woainal  System 

Proa  the  block  diagram  of  the  Pig.B.1  in  the 
appendix  5  one  can  obtain  the  following  nominal  state  equa¬ 
tions  of  the  uncoupled  pitch  autopilot. 


X1  »  C2  A3  X2  ♦  C2  A4  X3 


(3.1) 


X2 


XI 


K  Cl  A1  X2  -  K  Cl  A2  X3 


(3.2) 


X3  «  -  C3  X3  ♦  C3  Conv  X6 


(3.3) 


XU  «  -  Cl  C4  A1  X2  -  Cl  C4  A2  X3  -C4  X4 


(3.4) 
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X5 


C7  XU  -  C5  X5  -  C6  SZC 


(3.5) 


X6  »  C9/ConT  ♦  C2  C8/conT  A3  X2 


(3.6) 


♦  C2  C8/Conr  AU  X3  -  C7  C8  X4  ♦  (C5  C8  -  C9)  X5 

♦  C6  C8  MZC 


For  the  purpose  of  using  the  linear  nethod  these 
state  equations  oust  be  presented  in  matrix  form  as  given  by 


1  -  *<1^1  -  K 


O  -C3 


o  o  o 

ce  CtCq  A3  ctCe  a^ 

conv  eonv  conv 
V  - 


(3.7) 


The  correspondence  of  the  state  vectors  is: 

t 

XI  =  c|  ,  X2  »  ^  ,  13  «  Jp  ,  X4  *  X  ,  X5  =  I  ,  X6  = 

Definition  of  the  constants  Cl  through  C9  are  given 
in  appendix  B.  The  parameters  of  interest  for  this  system 
are  given  by  A1,  A2v  A3,  and  A4  which  are  : 


A1 


A2 


A3  «  ,  A4  *  . 


•••»«,  "'ap 

This  nomenclature  is  used  here  to  easily  apply  the 
sensitivity  theory  and  to  avoid  the  inconvenience  of 
carrying  symbols  and  constants. 
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2 .  se^si^iZltz  Eqaatians 

To  apply  the  procedure  developed  in  chapter  2, 
consider  the  trajectory  sensitivity  equation  (Eqa.2.13). 


\  Lu(t)  »  Aj^'‘o)=o  (3.8) 

■»  "  ■»  'b'Xj  oto° 


Froa  Bqn. 

3.7  one  can  see 

that 

o 

o 

o  ^ 

1 

-  H^C|  Af 

-KCfA^ 

o 

O  O 

o 

O 

o 

o  CjConv 

A  = 

o 

o 

o  o 

o 

o 

o 

Cf 

1 

O 

C9 

conv 

c,ce  As 
conv 

c»ce  A^ 

cowv 

*  ^Ss 

'sV'b  o 

k 


i>r  [  o 


O  o  O  -Cq 


(3.10) 


The  partial  derivatives  and  --g— -l^^^are 
evaluated  considering  the  paraaeters  of  interest  &1,  12,  A3 
and  A4,  that  are  respectively  the  aerodynaaic  coefficients 
present  in  the  pitch  channel. 
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Applying  the  partial  derivatives  with  respect  to  the 
paraeeters  we  have  the  following  aatrices 


o  o 

I 


o 

Conv 


(3.11) 


o  o 


"iA  Jo 


o  o 


(3.12) 


o  O  -KC,  O 
~  o  O  O  o 
'DA3  O  o  o 

00  o  o 


(3.  13) 


.  9  o  o 

o  o 

O  o  o 


(3.  14) 


and 


-251,..  .  0 

because  B  is  independent  of  the  paraseters  of  interest. 
In  teens  of  coaponents  one  can  see 


1- 


(3.15) 


i=1,2,..  .,6  j»1,2,3  and  4 

Here*  four  nodels  are  necessary  to  study  the  effect 
of  parameter  variatiens.  These  models  are  shown  in  Fig. 3.1. 

For  instance,  when  one  apply  the  theory  for  one 
parameter  of  interest,  say  &1,  the  sensitivity  equations  can 
be  obtained  from 


X 

m  m 

V  * « •  X, . 

!-  - 

r 

N 

VH 

= 

It 

•  0 

^  • 

• 

+ 

U^i«l 

•  <  > 

xs 

(3.16) 


Similar  procedure  can  be  done  for  the  other  parane* 
ters  12,  13,  14. 
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NoMiMAt  Hopei.  I  ^ 


C.  SEISITITITT  ABALZSIS 

A  coaputer  prograa  using  4]  vas  vrittan  in  order 

to  slaulate  slaultaneously  the  noeinal  and  the  sensitivity 
Bodels  as  shovn  in  the  appendix  0.  A  step  was  applied  as 
input  of  the  noeinal  systea. 

The  nuaber  of  eguations  solved  are  6  for  the  noainal 
aodel  and  24  for  the  sensitivity  aodel.  Each  piraaeter  was 
varied  slaultaneously  10X  froa  all  the  noainal  values. 


D.  AHALTSIS  OF  -PABABETEB  TABIAIIOHS 


The  results  are  plotted  in  Fig.  3.2  to  3.6  and  Table  I 
Each  state  variable  output  is  plotted  separately  with  the 
correspondent  four  output  sensitivity  functions.  By  aeans  of 
these  plots, the  following  observations  can  be  aade: 

The  plot  of  — U*  '  and  4  in  fig. 3.2 

indicates  that  a  paraaeter  change  a  Aj^(j=1,2, 3,  and  4) 
priaarily  affect  the  rise  tiae  and  overshoot  of  II  (t)  since 
4.11#  4.12,  4.13,  4.14,  are  largest  at  the  tine  where  these 

effects  in  XI  occur.  ^ 

The  plot  of  4.a^*  i*^»2,3,  and  4  in  Pig.3.3 

indicates  that  Z2(t)  is  little  affected  due  to  A  A|o  and 
^  A20  4nd  strongly  affected  in  the  rise  tine  due  to  paraa¬ 
eter  variations  4^  A30  and  A  a^o*  overshoot  and  steady 

state  are  little  affected  doe  to  paraaeter  variations 
(j-1#2,3,  and  4). 

The  plot  of  4.31*  ^*^'2,3,  and  4  in  Fig.3.4 

Indicates  that  I3(t)  is  little  affected  in  ^he  rise  tiae  due 
to  paraaeter  variations  A  (j«1,2,3,and  4).  The  overshoot 
and  steady  state  are  little  affected  due  to  paraaeter  varia¬ 
tions  4AA10  and^Ajo  and  strongly  affected  due  to 

AAj^. 
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The  plot  of  and  4  in  Pig.  3.5 
indicates  that  X4(t)  is  little  affected  in  the  rise  tine  and 
overshoot  dne  to  paraeeter  variations  ^&io  and  ^  A  20  and 
strongly  affected  dne  to  a  and  ^  A^. 

The  plot  of  and  4  in  Fig.  3.6 
indicates  that  X5(t)  is  little  affected  in  the  rise  tiee  due 
to  paraseter  variations  A  over* 
shoot  is  little  affected  due  to  A  A|^  and  A  k^^^  and  strongly 
affected  due  to  ^A^q  and  A.  A^^. 

The  plot  of  ^*^»2,3,  and  4  in  Fig. 3.7 
indicates  that  X6(t)  is  little  affected  in  the  rise  tine  due 
to  paraeeter  variations  A  Aj^  (^^1,2,3,  and  4).  The  over¬ 
shoot  and  steady  state  are  little  affected  due  to  A  A,^  and 
A  A20  strongly  affected  dne  to  A  a^^  and  ^  &40* 

Table  I  shovs  the  above  analysis  in  a  condensed  vay  to 
give  a  general  picture  of  all  states  and  sensitivity  func¬ 
tions  with  the  correspondent  effect  as  a  function  of  tine. 


B.  PABAHSTBB-IBDOCBD  OOTPOT  ABALISIS 

If  A«<  «  Xo  #the  actual  output  can  be  written  as 
(Bgn.2.29) . 


y(trX)  «  y(t,o^)  ♦  G"(t,o^)Aoc 


(3.17) 


The  coaputer  prograa  as  shown  in  the  Appendix  D  gives 
the  actual  output  when  10Xr30!(  and  40X  of  variation  froa  the 
noainal  value  of  each  paraeeter  is  assuaed.  Figs.  3.8 
through  3.13  show  the  actual  output  when  each  paraaeter  is 
varied  10%  froa  the  noainal  value. 

Fig. 3. 8  indicates  that  the  overshoot  and  rise  tiae  of 
the  actual  and  noainal  output  of  XI  are  strongly  affected 
which  is  in  agreeaent  with  the  previous  analysis. 
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Pig. 3. 9  iadicates  that  the  steady  state  and  overshoot  of 
the  actual  and  noeinal  output  of  Z2  are  priearlly  affected 
which  is  in  agreeaent  with  the  previous  analysis. 

Fig. 3. 10  indicates  that  the  staady  state  and  overshoot 
of  the  actual  and  nominal  output  of  Z3  ara  primarily 
affected  which  is  in  agreeaent  with  the  previous  analysis. 

Fig. 3. 11  indicates  that  the  overshoot  of  tha  actual  and 
nominal  output  of  Z4  is  primarily  affected  which  is  in 
agreeaent  with  the  previous  analysis. 

Fig. 3. 12  indicates  that  overshoot  of  the  actual  and 
nominal  output  of  Z5  are  primarily  affected  which  is  in 
agreeaent  with  the  previous  analysis. 

Fig. 3. 13  indicates  that  the  staady  state  and  overshoot 
of  the  actual  and  nominal  output  of  16  are  primarily 
affected  « having  little  effect  in  the  rise  time  which  is  in 
agreement  with  the  previous  analysis. 

Figs. 3. 14  through  3.19  and  Figs. 3. 20  through  3.25  show 
respectively  the  actual  output  for  SOU  and  40%  of  the 
nominal  value. 

From  the  plots  one  can  see  that  for  small  parameter 
variations  the  parameter-induced  output  error  is  negligible 
and  when  the  variation  becomes  large  as  30%  or  4  0%  one  note 
that  the  error  becomes  pronounced  and  that  modeling  is 
starting  to  break  down.  This  agrees  with  the  assumption 
made  in  the  derivation  of  the  Eqn.2.28. 

F.  COOPIBO  BOLL-IiB  iOTOPItOT  iliLZSZS 

1.  UAfiftc  ^hs  woi^ftaA  system 

From  the  block  diagram  of  the  Fig.B.1  in  the 
appendix  B  one  can  get  the  following  state  equations  of  the 
coupled  roll-yaw  autopilot. 
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Z1  «  B  conv  (&4  X3  ♦  i3  18  ♦  &5  Z11) 


(3.18) 


X2  »  C  Conr  (A7  X3  ♦  48  X8  ♦  A6  XII) 


(3.19) 


X3  «  -  XI  -  (4LPH4B/ConT)  X2  (3.20) 

♦  KB  4  42  X3  ♦  KB  4  41  XII 

X4  »  -  8  X4  -  17.6  X12  ♦  17.6  PHC  (3.21) 


X5  «  -  (0.755/CoaT)  X2  -  C  D  47  13  (3.22) 

♦  (0.755  -  8  D)  X4  -  5  X5  -  C  D  48  X8  -  C  D  46  XII 
-  17.6  X12  +17.6  PHC 

X6  «  C  E  47  X3  -  6  X6  ♦  B  C  48  X8  ♦  C  B  46  XII  (3.23) 


X7  «  -  F  KC  (0.755/ConT)  X2  (3.24) 

-(D  ♦  E)  F  KC  C  47  X3  ♦  F  KC  (0.755  -  8  D)  X4 
♦  KC(15  -  5  F)  X5  ♦  KC(6F  -15)  X6  -  15  X7 

-  (D  ♦  E)  F  KC  C  48  X8  -  (D  ♦  E)  F  KC  C  46  XII 

-  F  KC  D  17.6  X12  ♦  F  KC  D  17.6  PHC 

Z8  «  188.4  Coav  X7  -  188.4  X8  (3.25) 
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X9  «  (K1  A  A2/t1)  XI  -  X9  /  T1  ♦  (K1  A  A1/T1)  111 


XI 0  «-  (K2/Con7)  XI  -  (K2  H  ALPBAB/Conv  X2 

(k2/10)  (K1  A  A2/T1)  ♦  B  A4  -  H  ALPHAS  C  A7  X3 
*  (K2/10)  (B  A3  -  H  ALPHAS  C  A8  X8  4-  K2(1 

-  1/(10  T1))  X9  ♦  (K2/10)((K1  A  AH /T1  ♦  B  A5 

-  H  ALPHAS  C  A6)  Xl  1 

XII  =  188.4  Conv  X10  -  188.4  XII 


XI 2  *  X2  /  Conv 

Eqn.3. 30  glvQS  the  aatrlz  representation 
state  variables  above  aentioned. 
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(3.  26) 

(3.27) 


(3.28) 

(3.29) 
of  the 


(3.30) 
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The  correspondence  of  the  states  is: 
X1»r,X2»p,X3  =  fh  , 

X4  »  X  ,  X5  »  II  ,  X6  =  XI  , 

,  X8  =  Stt.  ,  X9  *  I  , 

X10  =  &Yc  »  111  “  Sy  ,  X12  *  4> 

The  definition  of  the  constants  Cl  throjgh  C25  are 

given  in  appendix  B.  The  paraseters  of  interest  for  this 

systes  are  given  by  &1  through  &8  which  are  : 

i1  «  Cy,  ,  A2  »  ,  A3  =  ,  A4  =  C„  ,  A5  *  C  , 

•'Aft  ngy 

A6  *  C.  ,  A7  =  C,^  ,  A8  =  Cfl  . 

<4y 

2 .  Sensitivity  Bgaations 

As  showed  in  the  previous  analysis  here  one  can 
apply  the  sane  procedure  using  the  TRAJECTOBI  SENSITIVITY 
EQUATION  as  given  in  Equ.2. 13. 


'Da 


L^o  ♦ 
Vo 


V, 


(3.31) 


Eor  the  purpose  of  this  procedure  again  one  can  see 
in  Eqn.3.30  the  correspondent  ■atri£es  A^  and  ^ 

The  partial  derivatives 

evaluated  considering  the  paraseter  A  of  interest.  In  this 
case  they  are  A1 ,A2# ...» A8,  that  are  respectively  the  aero¬ 
dynes  ic  coefficients  present  in  the  coupled  roL 1-yaw  auto¬ 
pilot. 

Once  again r applying  the  partial  derivatives  with 
respect  to  the  paraseters  one  can  obtain  eight  satrices 
respectively  as  found  ,sinilarly»  in  the  previous  case  of 
the  uncoupled  pitch  autopilot.  ^ 

Bere»one  can  see  that  «  0  because  B  is  inde- 

pendent  of  the  paraseters  of  Interest. 

In  terns  of  cospooents  one  obtains 
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6.  SBVSITITZIT  MALTSIS 

In  ord«r  to  slaolate  siaoltannoaly  the  noaiaal  and  the 
sensitivity  sodels^a  cospnter  prograa  was  written  as  shown 
in  the  appendix  B. 

For  analysis  par pose, each  paraseter  was  varied  slaalta- 
neously  10X  from  the  noaiaal  valae. 


-f 
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1.  InalTSis  of  oc-Pac^aetes 

X 

The  plots  of  A,;  ®  -^1“*: — Ia'  »  j*1,2,...,8  in 

Pig.  3.27  and  3.20  indicate  that  a  parameter  change 
(j*1 »  2,.  ..  ,  8)  primarily  affect  the  overshoot  of  XI  (t). 

The  plots  of  A^j‘  =  j»1,2,...,8  in 

Fig.  3.29  and  3.30  indicate  that  a  parameter  change 
(j*1 »  2,. .. ,  8)  primarily  affect  the  overshoot  of  X2(t). 

The  plots  of  Ajj-  s  j=1,2,...,8  in 

Pig. 3.31  and  3.32  indicate  that  a  parameter  change  AAJq, 
<js1 , 2,. .. ,  8)  primarily  affect  the  overshoot  of  X3(t). 

The  plots  of  A^-  a  irA^“”*A3o'  j=1f2,...,8  in 

Fig. 3.33  and  3.34  indicate  that  a  parameter  change 
(j*1  #  2 ,  8)  primarily  affect  the  overshoot  of  X4{t). 

The  plots  of  Agj  *  j=1.2,...,8  in 

Pig. 3.35  and  3.36  indicate  that  a  parameter  change  AAj^, 
(jal  »2,...  ,8)  primarily  affect  ^e  overshoot  of  X5(t). 

The  plots  of  j=1,2,...,8  in 

Pig.  3.37  and  3.38  indicate  that  a  parameter  change 
(jal  ,2/. ..  r8)  primarily  affect  the  overshoot  of  X6  (t) . 

The  plots  of  Ayj*  =  »  j=1»2,...,8  in 

Pig.  3.39  and  3.40  indicate  that  a  parameter  change 
(j*1 » 2,.  . . ,  8)  primarily  affect  the  overshoot  of  X7(t). 

The  plots  of  Asii  “  6Aj””*A3o'  j=1»2,...,8  in 

Pig.  3.41  and  3.42  indicate  that  a  parameter  change  AAj'^, 
(j*1  #2,. ..  ,8)  primarily  affect  the  overshoot  of  X8(t>. 

The  plots  of  Agj  ®  j*1»2,...,8  in 

Pig.  3.43  and  3,44  indicate  that  a  parameter  change  AAJq, 
(j*!  »2,. ..  ,8)  primarily  affect  the  overshoot  and  rise  time 
of  19  (t)  . 

i  Aj  O 

The  plots  of  Aioj’  •  j»1,2,...,8  in 

Pig. 3.45  and  3.46  indicate  that  a  parameter  change 
(j*1  r  2». ..  r  8)  primarily  affect  the  rise  time  of  X10(t). 
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The  plots  of  -^1 


^  X  (I 

j*1r2,...,8  in 

Fig.  3.47  and  3,48  indicate  that  a  parameter  change  A3  o, 
» 2  ,  8)  primarily  affect  ^he  overshoot  of  X11(t). 

The  plots  of  A 13^'  =  *-^00'  j=1»2,...,8  in 

Pig.  3.49  and  3.50  indicate  that  a  parameter  change 
(j*^ » 2  ,. . . ,  8)  primarily  affect  the  overshoot  of  Xl2(t). 

Table  II  and  III  show  the  above  analysis  in  a 
condensed  way  to  give  a  general  picture  of  all  states  and 
output  sensitivity  functions  with  the  correspondent  effect 
as  function  of  time. 


2.  Parana  ter- Induced  Output  Analysis 

As  Shown  previously,  if  Aoc  «  ot.^  , 
output  can  be  written  as 


the  actual 


(3.34) 


The  computer  program  given  in  appendix  B  was  written 
for  simulating  the  system  when  10%,  30%,  and  4D%  of  varia¬ 
tion  from  the  nominal  value  of  each  parameter  is 
assumed. Pig. 3. 51  through  3.62  give  the  actual  output  when 
each  parameter  is  varied  10%  from  the  nominal  value. 

Pig. 3. 63  and  3.64  give  the  actual  output  for  30%  of 
variation  from  the  ncninal  value  of  each  parameter  assumed. 
Pig. 3.65  and  3.66  give  the  actual  output  for  40%  of  varia¬ 
tion  from  the  nominal  value  of  each  parameter  assumed. 
These  plots  show  the  output  of  the  state  variables  13  and 
111  that  present  strong  deviations  just  to  give  the  behavior 
of  the  system  when  parameter  variations  are  not  small.  One 
notes  that  modeling  is  starting  to  break  down. 
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figar«  3.5  S«n«ltiTity  of  m  with  Respect  to 
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Figar*  3.8  ictoal  and  loainal  oatpot  of  xi  (lOS  Tarlation) 


rlgar*  3.10  ictaal  and  Roalnal  Ootpat  of  Z3  (lOX  ▼ariatlon) 


Pigur*  3.13  Actual  and  Soninal  Output  of  X6  (13X  Tariation) . 
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Flgor«  3.15  ictaal  and  Honinal  Oatpat  of  12  (33%  variation) 


Plgur«  3.16  Actaal  and  Hoalnal  Oatpat  of  X3  (30X  variation) . 
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Plgare  3.17  ictaal  and  Noainal  Output  of  X4  (3DX  Yariation) 


W.gar»  3.I8 


Actual  aad  loaiaal 


ot  IS  „»*  _ 
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Plgare  3.19  Actual  and  Noainal  Output  of  Z6  (39%  variation)  . 
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Figure  3-20  ictual  and  foninal  ootpat  of  XI  («3f  variation). 


Figar*  3.21  ictaal  and  Hoalnal  oatpat  of  X2  (a9%  variation) 


rigor*  3.23  ictoal  and  Roainal  Ootpot  of  Z4  (HOX  Tarlatlon) 
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Flgar*  3.25  ictaal  and  Soalnal  Oatpat  of  X6  («3X  ▼ariatlon) 


Pigtir*  3.26  Boll-Ta«  Moainal  aad  Sensitivity  Bodsls. 
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Figur*  3.27  S^nsitiTity  of  I1  with  Raspect  to  k^,k2,k3,k^^ 


Pigar*  3.29  Sensitivitj  of  X2  with  Hespect  to  K^,k2,K3,M 


Plgare  3.30  Sensitivity  of  X2  with  Respect  to  i5,i6,i7,A8 


Flgar*  3.38  SensltiTlty  af  X6  wltti 
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Pignr*  3«44  S^asltiTltj  of  Z9  with  Bospoct  to  &5«A6,A7, i8 
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Figvr*  3.«S  S^nsltivltr  of  Z10  with  Bospoct  to 
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Pigar*  3.51  ictaal  and  Roalnal  Output  of  xi (10X  rarlation) . 
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Plgar*  3.52  Ictaal  and  Noainal  Oatpat  of  Z2(10X  Tarlatlon)  . 
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'O.OO  0.40  0.80  1.20  1.60 


Pigac*  3.55  Actaal  and  Moainal  oatpat  of  Z5(10X  Tariation) 


Figar*  3.58  Ictaal  and  Honinal  Oatpot  of  Z8  (10X  variation) 
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Tlgar*  3.62  Actual  and  loainal  output  of  XI 2 (10 X  tariatlon) 
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Pigar«  3.63  ictaal  and  lonlnal  Oatpat  of  X3  (301  Tarlatlon) 


Figar*  3.64  ictaal  and  lonlnal  Oatpat  of  Z11(30S  ▼ariatlon) 


Fig\ir«  3.65  ictaal  aad  Hoaiaal  Oatpat  of  Z3  (%0t  Taziation) 


riqnf  3.66  ictaal  and  Soninal  Oatpnt  of  X11{43S  variation). 


IV-  iPPLICATIOH  OT  SBBSITITITT  MAtlSIS  TO  N0HiIgE4R  ST  ST  BBS 


A.  IHTBODOCTIOH 

The  three  diaenslonal  nonlinear  system  used  here  for 
sensitlrlty  analysis  purposes  is  presented  in  appendix  C. 
The  eisslla  configuration,  size  and  mass  properties  are  the 
same  as  those  used  in  the  linear  system  and  were  presented 
in  appendix  &. 

In  section  B  the  sensitivity  agnations  with  respect  to 
parameters  of  nominal  system  is  shown  and  in  section  C  the 
sensitivity  equations  with  respect  to  parameters  that  showed 
to  be  of  more  effect  in  the  previous  analysis  is  presented. 

Section  D  shows  the  results  of  the  trajectory  sensi¬ 
tivity  equations  when  step  inputs  are  applied  at  specific 
time  as  mentioned  in  appendix  C. 

Section  E  gives  the  parameter-induced  outpjt  analysis 
when  each  parameter  is  varied  lOf  from  the  nominal  value. 

B.  lOILIBBAB  BQOATIOIS  OF  THE  NOHINAL  STSTBH 

Prom  the  bloch  diagram  of  Fig.  C.1  and  Bqn.  C.3  through 

C. 21  as  given  in  appendix  C  one  have  the  following  nominal 
equations. 


XI  *  -  C4  XI  -  C02  C4  K2 


(4.1) 


X2 


C5  X2  -  C6  HZC  ♦  C6  C19  C03(X7)  ♦  C7  XI  (4.2) 
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13  »  (C5  C8  -  C9)  X2  ♦  C6  C8  MZC  (4.3) 

-  C6  C8  C19  Cos(X7)  -  C7  C8  XI  ♦  (C8/Conv*)  XI 7  X18 

♦  Cl  C8  A1  ♦  (C9/ConT)  X5 

« 

X4  =  -  C3  X4  ♦  C3  Conv  X3  (4.4) 

« 

X5  «  (X17  X18)/Conv  ♦  Cl  Conv  A1  (4.5) 

X6  =  X5  -  KB  C02  A2  -  (X16  X18)/C3Q\r  (4.6) 

X7  *  X5  Cos(X19) /Conv  -  X17  Sin  (XI  9) /Conv  (4.7) 

X8  =  -  CIO  X8  ♦  (C12  -  C11  C13)  XIO  (4.8) 

♦  C11  C14  PHC  -  C11  C14  X19  -  C11  CO  (A8  X16  ♦  A7  X15 

♦  A3  X13)  -  (C12/ConT)  X18 

X9  *  -  X9/T1  ♦  K1/T1  C02(AA  X16  ♦  A6  X15)  (4.9) 

XIO  =  -  C13  XIO  ♦  C14  PHC  -  C14  X19  (4.10) 

XII  »  -  C15  XII  ♦  C16  C0(A8  X16  ♦  A7  X15  ♦  A3  X13)  (4.11) 
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X12  =  -  C17  X12  ♦  K(C17  -  (C10/C18)  X8 


(«. 12) 


(K/C18)  (C12  -  C11  C13)  X10  ♦  K  C11  (C14/C18)  P  HC 
-  K  C11  (C14/C18)  X19  -  (K/C18)  30  A8<C11  -  (K/C18) 

CO  A3(C11  ♦  C16)  X13  -  (K/C12)  (Cl  8  Conv)  X18 

♦  K(C15/C18)  -  C17)  XII 

X13  =  -  C3  X13  ♦  C3  conv  X12  (4.13) 


X14  *  (K2/10)  (-  (X9/T1)  C02(AA  X16  ♦  A6  X15)  (4.14) 

-  X15  (Xia/Conv*)  ♦  C01(A9  X16  ♦  A4  X15  ♦  A5  K13) 

-  (KTP/Conv*)  (X5  -  KB  C02  A2  -  XI 6  (X 1 8/Con v)  X18 

-  (KTP/Conv)  XX6  C0(a8  X16  ♦  A7  X15  ♦  A3  X13) 

♦  (K2/Conv)  X17  -  K2(KIP/ConT*)  XX6  X18 


X15  «  -  C3  X15  ♦  C3  Conv  X14 


(4.  15) 


X16  =  KB  C02(AA  X16  ♦  A6  115)  ♦  X5(X18/Conv)  -  X17  (4.  16) 


X17  =  -  (X5  XI 8) /Conv  ♦  C01  Conv(A9  X16  (4.17) 

♦  A4  X15  ♦  A5  X13) 

X18  a  CO  Coay(A8  X16  ♦  A7  X15  A3  X13)  (4.18) 


X19 


X18/Con7 


(4.  19) 
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The  correspondence  of  nonlinear  state  sectors^  as  given 
in  the  above  equations  with  relation  to  the  nonlinear  system 
presented  in  Fig.C. 1  in  appendix  C  are  : 

XI  =  X  ,  X2  »  I  ,  X3  =  ,  xi*  =  Sp  ,  X5  =  q  , 

X6  =  ,  X7  »  6  ,  X8  =  T1  ,  X9  =  12  ,  XI 0  *  X 1  , 

XII  =  12  ,  112  »  Sue,  X13  =  Sr  ,  X14  *  6y<.,  X15  =  6y  , 

X16  =  (i  ,  X17  »  r  ,  X18  *  p  ,  and  X19  *  4  • 

The  paraneters  of  interest  for  this  nonlinear  systen  are 
given  by 

41  A2  A3  AH  r  and  A5 

Vo. 

Definition  of  the  constants  CO  through  C19  are  given  in 
appendix  C. 

4s  previously,  this  noeenclature  is  used  here  to  easily 
apply  the  sensitivity  theory. 

C.  S0HLIHB4B  SEBSITITITT  BQ04TI0HS 

Proa  the  sensitivity  theory  given  in  chapter  2  one  knows 
froa  Bqn.2.12  that  for  nonlinear  systems  one  has 

^  'Dx® 

Here 

X  s  f  (X1,X2,.  .  .,X19,t,NZCrA1,42r..  .  rA5} 


^The  state  vectors  and  paraneters  in  the  nonlinear  case 
have  no  correspondence  to  the  state  vectors  and  parameters 
used  in  the  linear  case. The  correspondence  can  be  found  la 
the  appendices  B  and  C. 


(4.20) 


(4.21) 


Froa  Eqa.2.14  one  knows  that  tha  trajactory  ssnsitiTlty 
vector  is  given  as 


•"N  /  __ 

- - 1  .  A  (0)  «  0 

0  0^  oT® 


(4.22) 


Applying  the  above  theory  one  obtains  tha  nonlinear 
sensitivity  equations  as  given  in  Eqn.4.24  throagh  &4e42, 
where  W4S  considered  the  paraaetar  to  be  varied. 


ail  «  -  CUA11  -  C02  C4(DA26a.61  *■  DA24a41) 


(4 . 23) 


i.21  »  -  C5a.21  -  C6  C19  Sin(l7)  ^71  ♦  C7^11  (4.24) 

iai  =  (C5  C8  -  C9)  A21  ♦  C6  C8  C19  Sin  (X7)  ^71  (4.25) 

-  C7  C8  a.11  ♦  C8(x18A171  ♦  (X17  a.181)/(Conv*» 

♦  Cl  C8(DA16/i61  ♦  DA14  .141  ♦  1.)  ♦  (C9  .15 1) /Conv 

.141  »  -  C3/141  ♦  C3  Conv  yl31  (4.26) 

a.51  »  (1174.181  ♦  118^117 1)/Conv  (4.27) 

♦  Cl  Conv(DA16  4.61  ♦  DA14/141  ♦  1.) 

/161  *4.51  -  KB  C02(DA26/i61  *  0A24  4.4  1)  (4.28) 

-  ((X16/1181  ♦  X18  4161) /Conv 
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-  I 


i71  *  (Cos (XI 9)  >151  -  X5  Sic  (X19>  ^191) /COOT  (4.29) 

-  (Sin(Xl9)  ^171  ♦  X17  Cos(X19)  >i171)/ConT 

jl81  *  -  CIO  >181  ♦  (C12  -  C11  C13)  /1101  (4.30) 

-  C11  C14A191  -  C11  CO  (18  >1161  ♦  17>1151  ♦  k3^131 

♦  /161(X16  Da8  ♦  115  Da7  ♦  X13  D113))  -  (Cl  2  >11 81) /Conv 

>191  «  ->191/T1  ♦  K1/(T1  C02)  (11x1161  (4.31) 

*  16A151  *■  a61(X16  Dll  *  X15  Da6)  ) 

>1101  =  -  C13>1101  -  C14>1191  (4.32) 

/ 


>1111  *  -  C15/1111  ♦  C16  C0(18>1161  ♦  17  >1151  (4.33) 

♦  13  >1131  ♦  >161  (X16  D18  ♦  X15  D17  ♦  X13  D13)  ) 

>1121  =  -  Cl7xi121  ♦  K(C17  -  C10/C18)  >181  (4.34) 

♦  (K/C18)  (C12  -  C11  C13)  >1101  -  K  C11  (C14/C18)  xl191 

-  (K/C18)  C0(C11  ♦  C16)  (18/1161  ♦  17  >1151  ♦  13.1131 

♦  13  >1131  ♦/a61(X16  D18  ♦  X15  D17  ♦  X13  Da3)  ) 

-  K  C12/(C18  ConT)>1181  ♦  K  (C15/C18  -  C17)>1111 

• 

>1131  =  -  C3>1131  ♦  C3  Coqt/1121  (4.35) 
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(4.36) 


il41  »  (K2/10)  (-  X9/T1  ♦  (K1/T1)  Z02{kk 

a.161  *  &6/1151  *  ^61(X16  D&&  ♦  X15  0X6))  -  (X15/1181 

♦  X18  i151)/(CDn7*)  ♦  C01(A9A161  ♦  A4  2.151  ♦  A5  2.131 

♦  ^61(X16  DA9  ♦  X15  DA4  ♦  X13  DA5>  )  -  KYP/(C0!1V*) 

(15  2.181  ♦  X182.51  -  KB  C02(A2  2.181  *  Xl8(DA2Syd61 

♦  DA24/i41)  -(X18*ai61  ♦  2  X18  X162.181)  /CoaT)) 

-  (KlP/Conv)  C0(  zx6  2.161  ♦  X16  2.161)  ♦  A7  2.161) 

♦  A7(XX6  2.151  ♦  X152.161)  ♦  A3(XX6/i131  ♦  X13  2.161) 

♦  X16(XX6  DA8  2.61  ♦  A8  2-161)  ♦  X15(XX6  DA7/161 

♦  A7  2161)  *  X13(XX6  DA3  261  +  A3  2161))) 

♦  (K2/Contr)  2171  ♦  K2  291  -  K2  KrP/(Conv*) 

(XX6  2181  *  X18  2161) ) 

2151  =  -  C3  2151  ♦  C3  CoaT2141  (4.37) 


2161  »  KB  C02(AA  261  ♦  A6  2151  ♦261(  («».  38) 

X16  OAA  ♦  X15  DA6))  ♦  (X62181  ♦  X18  261)/C0Qf  -2171 
2171  =  -  (X5  2181  ♦  X18251) /Conv  (4.39) 

♦  C01  ConY(A9  2161  ♦  A4  2151  ♦  A5  2131  ♦261(116  DA9 

♦  XI 5  0A4  ♦  X13  DA5)) 

2181  =  CO  cony  (A8  2161  ♦  A7  2151  (4.40) 

♦  A3  2131  ♦261(X16  DA8  ♦  X15  DA7  ♦  X13  DA3)) 


2191  »  21 81/Cony 


(4.41) 


The  saee  procedure  abotre  east  be  eade  for  the  other 
paraeetars  that  eost  affected  the  pratious  analysis  and  will 
be  giwen  in  the  nonlinear  sensitiwity  analysis  presented  in 
the  next  section. 


0.  lOILIlBUt  SBSSITITITT  iBiLTSIS 


k  coeputer  prograe  in  appendix  3  shows  the  simulation  of 
the  sensitiwity  equations  with  respect  to  the  selected 
paraaeters. 

The  nueber  of  equations  solved  are  19  for  the  nominal 
equations  and  95  for  the  sensitivity  equations.  Here,  the 
paraaeters  were  selected  from  the  previous  two  cases  given 
in  Chapter  III  that  showed  to  have  most  effect  in  the  time 
response  of  the  systems. 

Proa  the  analysis  performed,  two  paraaeters  were 
selected  from  the  uncoupled  pitch  autopilot  and  three  from 
the  coupled  roll«yaw  autopilot,  which  has  been  shown  to  be 
more  sensible  to  parameter  variations.  These  paraaeters  are 
respectively  (A2)  and  0,^^^  (A4)  in  the  'incoupled  pitch 
autopilot  andQj^^(A3),  ,  and  C,,^^(A8)  in  the 
coupled  roll-yaw  autopilot. 

In  order  to  compare  the  influence  of  the  parameter  of 
interest  in  the  nonlinear  case  some  state  variables  were 
chosen  to  be  analysed.  The  state  variables  selected  were: 


From  the  uncoupled  pitch  autopilot  : 

X3  »  ,  X4  «  Sp  ,  X5  »  q  ,  X6  » 

From  the  coupled  roll-yaw  autopilot  : 

X12  =  ,  X13  -  Sr.  ,  X14  -  ,  X15  -  Sy  , 

X16  »  ,  X17  »  r  ,  X18  «  p  ,  and  X19  »  ^  . 
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The  results  are  plotted  in  Fig.U.  1  through t  4.12  and 
Table  17  and  7.  Each  state  variable  output  is  plotted 
separately  with  the  correspondent  five  output  sensitivity 
functions. 

From  the  plots,  the  following  observations  can  be  sade: 

Fig. 4.1  indicates  that  the  state  variable  is 

little  affected  in  the  rise  tise  by  the  paraosters  hi  and 
h2.  The  overshoot  is  strongly  affected  by  A1  and  &4  and 
little  affected  by  A2  and  A3.  The  steady  state  is  strongly 
affected  by  A1  and  little  affected  by  A2. 

Fig. 4. 2  indicates  that  the  state  variable  X4(£p)  is 

little  affected  in  the  rise  tine  by  A1  and  A2.  The  over¬ 
shoot  is  strongly  affected  by  A1,  A3,  and  A4  with  little 

effect  by  A2  and  AS.  The  steady  state  is  srrongly  affected 
by  A1  and  little  affected  by  A2. 

Fig.  4. 3  indicates  that  the  state  variable  X5  (  )  is 

strongly  affected  in  the  rise  tine  by  A1  and  A2.  The  over¬ 
shoot  is  strongly  affected  by  A3,  A4,  AS,  and  little 

affected  by  A1  and  A 2.  The  steady  state  is  not  affected  by 

paraaeter  variations. 

Fig.  4. 4  indicates  that  the  state  variable  X6  (  «<  )  is 

strongly  affected  in  the  rise  time  by  A1  and  A2.  The  over¬ 
shoot  is  strongly  affected  by  A3,  A4,  and  AS.  The  steady 

state  is  little  affected  by  A1  and  A2. 

Fig.4.S  indicates  that  the  rise  tiae  and  steady  state  of 
the  state  variable  X12(if^  are  not  affected  by  paraaeter 
variations.  The  overshoot  is  strongly  affected  by  A3  and  AS 
and  little  affected  by  A1,  A2,  and  A4. 

Fig. 4. 6  indicates  that  the  rise  tiae  and  steady  state  of 
the  state  variable  X13(^K^)  are  nst  affected  by  paraaeter 
variations.  The  overshoot  is  strongly  affected  by  A3,  A4, 

and  AS  and  little  affected  by  A1  and  A2. 

Fig.  4. 7  indicates  that  the  rise  tiae  and  the  steady 
state  of  the  state  variable  114  not  affscted  by  the 
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pacaaetar  Tarlatioas.  The  07ersfaoot  is  strongly  affected  by 
&4  and  &5  and  little  affected  by  \2,  and  &3. 

Pig. 4. 8  indicates  that  the  rise  tine  and  the  steady 
state  of  the  state  variable  X15(6y'  t  are  not  affacted  by  the 
paraaeter  variations.  The  overshoot  is  strongly  affected  by 
&3y  &4,  and  &5  and  little  affected  by  &1,  and  \2. 

Fig. 4. 9  indicates  that  the  rise  time  and  the  steady 
state  of  the  state  variable  Z16(^  )  are  not  affacted  by  the 
paraaeter  variations.  The  overshoot  is  strongly  affected  by 
A3,  A4,  and  AS  and  little  affected  by  A1,  and  A2. 

Fig. 4. 10  indicates  that  the  rise  tiae  and  the  steady 
state  of  the  state  variable  X17(  X"  )  QOi  affected  by  the 
paraaeter  variations.  The  overshoot  is  strongly  affected  by 
A3,  A4,  and  AS  and  little  affected  by  A1,  and  A2. 

Fig.  4. 11  indicates  that  the  rise  tine  and  the  steady 
state  of  the  state  variable  X18(  )  are  not  affected  by  the 
paraaeter  variations.  The  overshoot  is  strongly  affected  by 
A3,  A4,  and  AS  and  little  affected  by  A1,  and  A2. 

Fig. 4. 12  indicates  that  the  rise  tiae  and  the  steady 
state  of  the  state  variable  X19(  ^  )  are  not  affected  by  the 
paraaeter  variations.  The  overshoot  is  strongly  affected  by 
A3,  A4,  and  AS  and  little  affected  by  A1,  and  A2. 

E.  PABABETEB-IHOOCBD  OOTPOT  AHALTSIS 

As  shown  previousy,  if  ,  the  actual  output  can 
be  written  as 


y  (t,OC)  t  y  (t,  )  ♦  <r  {tf  Ot^)  £kO<. 


(4.42) 


The  coaputer  prograa  given  in  appendix  6  was  written  for 
slaulating  the  systea  when  variations  froa  the  noainal  value 


0 


/ 


3 


of  each  paraaeter  is  assuaed.  Figs. 4. 13  throaght  4.24  give 
the  plots  of  the  actaal  and  noainal  output  of  the  state 
variables  that  are  present  in  the  unclouped  pitch  autopilot 
and  coupled  roll-yaw  autopilot  when  10X  of  paraaeter  varia¬ 
tion  is  assuaed..  The  others  state  variables  ace  not  given 
here  because  they  are  just  interaediate  state  variables. 

these  results  conflra  the  theory  presented  in  chapter  2 
when  saall  paraaeter  variations  are  assumed.  Pigs. 4. 25  and 
4.26  present  the  actual  and  noainal  output  of  the  state 
variables  Z3  and  Z11  that  showed  pronounced  variations  when 
30%  of  paraaeter  variation  is  assuaed.  Figs.4.27  and  4.28 
present  the  actual  and  noainal  output  of  X3  and  Z11  when  40% 
of  paraaeter  variation  is  assuaed.  From  these  polts  one 
note  that  the  aodeling  is  starting  to  break  down. 
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Flgnc*  4*2  SvnsitlTlty  of  Z4  with  respect  to 
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Figure  4.4  SensitlTlty  of  16  with  respect  to  &1 ,i2,A3,A4,A5. 
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Pigar*  4.6  Sansititrlty  of  X13  vitli  raspact  to 
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figure  4.10  Sensitivity  of  Z17  with  respect  to  11,A2«13» 
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Figar*  4.11  sansitlTlty  of  X10 
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TIME 


Pigar«  4.14  Ictaal  and  Voainal  Output  of  X4  (18%  Tariation) . 
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Fifar*  4.15  Actaal  and  loainal  Oatpat  o£  X5  (10 S  Tariatlon) 
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Flgar*  4.17  tctaal  and  Soninal  Oatpat  of  112  (10X  variation). 
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Plgare  4.20  ictaal  and  Roainal  Oatpat  of  X15  (10%  Tariation). 


Pigat*  *-21  ictaal  aad  loainal  Oatpat  of  X16  (10*  rariatlon) 


ictttal  and  loniaal  output  of  X18  (10*  Tariation) . 
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Figure  4.25  ictual  and  Boalnal  Output  of  X6  (30S  Tariation) . 
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Influence  of  Paraaetecs  in  the  Piee  Response  (cont.) 
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?,  COHCIOSIOMS  MD  RBC3HIIBNDATI0HS 
A.  COSCLOSIOHS 

The  results  of  the  nonlinear  3-D  parameter  sensitivity 
analysis  presented  in  Chapter  IV  verify  the  linear  analysis 
of  Chapter  III.  These  conclusions  can  be  quickly  made, 
using  Tables  I,  II,  III  and  IV  that  give  a  brief  review  of 
the  time  response  of  the  linear  and  nonlinear  system. 
Another  means  of  comparison  of  the  linear  and  nonlinear 
analysis  is  to  use  the  figures  that  give  a  precise  visuali¬ 
zation  of  the  time  response. 

The  comparison  of  the  linear  and  nonlinear  analysis  can 
be  done  ty  means  of  the  plots  as  follows: 


B.  EFFECT  OH  STATE  VABIABLES  OOS  TO  PABAHETER  VARIATIONS 

1.  Effect  on  State  Variable  Due  to  Variation  in 

c  (i,i£)  : 

Figs. 3. 7  and  4.1  (  Ac'^  and  A^l  )  show  that  is 

little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot  and  steady  state. 

Figs.  3. 4  and  4.2  (  and  X4  j  )  show  that  is 

little  affected  in  the  rise  time  and  strongly  affected  in 
the  overshoot  and  steady  state. 

Figs. 3. 2  and  4.3  (  and  A/Ji)  show  that  g  is 

strongly  affected  in  the  rise  time,  little  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Figs. 3. 3  and  4.4  (  and  Agi  )  show  that  <<  is 

strongly  affected  in  the  rise  time,  little  affected  in  the 
overshoot  and  steady  state. 
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2»  Effect  on  State  Variable  Dae  ^  variation  in 

C_(Oirit)  : 

M 

Figs. 3. 7  and  4.1  (  and  X32)  show  that  ^Pc  is 

little  affected  in  the  rise  time,  overshoot  and  steady 

stat  e. 

Pigs.  3.4  and  4.2  (;^g2.  )  show  that  Sp  is 

little  affected  in  the  rise  time,  overshoot  and  steady 

state. 

Figs. 3. 2  and  4.3  (  A12  and  Aaz)  show  that  g  is 

strongly  affected  in  the  rise  time  and  little  affected  in 
the  overshoot  and  steady  state. 

Figs. 3. 3  and  4.4  (Azz  ^62)  show  that  is 

strongly  affected  in  the  rise  time  and  little  affected  in 
the  overshoot  and  steady  state. 

3.  Effect  on  State  Variable  Due  to  Variation  in  C . 

Figs. 3. 40  and  4.6  (A7g  and  XiZ3)  show  that  is 

little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Figs. 3. 42  and  4.7  (Xggand  A155)  show  that  is  not 
affected  in  the  rise  time  and  steady  state,  and  strongly 
affected  in  the  overshoot. 

Figs. 3. 46  and  4.8  (  AioQ  )  show  that  6y<^  is 

little  affected  in  the  rise  time,  little  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Figs. 3. 48  and  4.9  (  A|ig  and  X153  )  show  that  is 

little  effected  in  the  rise  time,  strongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Figs. 3. 32  and  4.10  (XsftandAi^^  )  show  that  is 
little  affected  in  the  rise  time,  stronglu  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Pigs. 3. 28  and  4.11  (^j^and  A\73 )  show  that  \r  is 
little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 


Figs. 3. 30  and  4.12  Al83  )  show  that  p  is 

littls  affectad  in  the  rise  time,  strongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Figs. 3. 50  and  4.13  (A(2g  show  that  ^  is 

not  affected  in  the  rise  time  and  steady  state  and  strongly 
affected  in  the  overshoot. 

4.  Effect  oii  State  yariable  Due  ^  Varia^on  in  C  _  ; 

Figs. 3. 40  and  4.6  and  Aj2.q)  show  that  is 

not  affected  in  the  rise  time  and  steady  state  and  little 

affected  in  the  overshoot. 

Figs. 3. 42  and  4.7  and  A154  )  show  that  6ix  is 

not  affected  in  the  rise  time  and  steady  state  and  strongly 

affected  in  the  overshoot. 

Figs. 3. 46  and  4.8  (A|q5-  and  A|^<^  )  show  that  is 
little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Figs. 3. 48  and  4.9  (Aji^  and  Ai5^)  show  that  Sy  is 
not  affected  in  the  rise  time  and  steady  state  and  strongly 
affected  in  the  overshoot. 

Figs. 3. 32  and  4.10  (As^and  A  14»4  )  show  that  p)  is 
little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Figs. 3. 28  and  4.11  (  Ais  and  Ajy^)  show  taat  r  is  not 

affected  in  the  rise  time  and  steady  state  and  strongly 
affected  in  the  overshoot. 

Figs. 3. 30  and  4.12  (AzS  and  Aig^  )  show  that  p  is 
little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Pigs. 3. 50  and  4.13  (A|2.5  A|9^)  show  that  ^  is 

not  affected  in  the  rise  time  and  steady  state  and  strongly 
affected  in  the  overshoot. 
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5»  Effect  on  State  Variable  Djig  to  Yariatioa  in  C  _  : 

“  5  ^ 

Pigs. 3. 39  and  4.6  (^Y^and  A|Z6‘)  show  that  Sf^^. 
little  affected  in  the  rise  tiae,  strongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Pigs.3.41  and  4.7  (^^3  and  A(35  )  show  that  is 

little  affected  in  the  risr  tiae,  stronglu  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Figs. 3. 45  and  4.8  and  )  show  that  6yc  is 

little  affected  in  the  rise  tiae,  strongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Pigs. 3. 47  and  4.9  (Ajia  ^ISS^  show  that  Sy  is 
little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Figs. 3. 31  and  4.10  (A33aai  AifcS  )  show  that  ^  is 
little  affected  in  the  rise  time,  stongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Pigs. 3. 27  and  4.11  (Ajj  and  A 175  )  show  that  }r  is 
little  affected  in  the  rise  tiae,  strongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Figs. 3. 29  and  4.12  (Azj  show  that  Ip  is 

little  affected  in  the  rise  tiae,  strongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Pigs. 3. 49  and  4.13  {Ai25  aoid  Aig^)  show  that  ^  is 

little  affected  in  the  rise  tiae,  strongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

C.  EFFECT  OF  TBUECTORI  SBHSITITITT  POHCTl-'HS 
1 .  0ncoup3,ed  Pitch  Aat opilot 

The  trajectory  sensitivity  functions  All,  /112, 
A13,  and  i-n  the  linear  case  and  the  correspondent 

trajectory  sensitivity  functions  151  and  152  ii  the  nonli¬ 
near  case  show  strong  effect  on  the  rise  tiae  aad  overshoot 
of  the  state  variable  q. 
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The  trajectory  sensitivity  functions  /12  3  and  /i24 
in  the  linear  case  and  the  correspandent  trajectory  sensi¬ 
tivity  function  161  in  the  nonlinear  case  show  strong 

effect  on  the  rise  time  of  the  state  variable  . 

The  trajectory  sensitivity  functions  /13  3  and  /i34 
in  the  linear  case  and  the  correspondent  trajectory  sensi¬ 
tivity  function  yi4  1  in  the  nonlinear  case  show  strong 

effect  on  the  overshoot  and  steady  state  of  the  state  vari¬ 
able  . 

The  trajectory  sensitivity  functions  A5 3  and  ^64 
in  the  linear  case  and  the  correspondent  sensitivity  trajec¬ 
tory  function  A31  in  the  nonlinear  case  show  strong  effect 
on  the  overshoot  and  steady  state  of  the  state  variable  . 

2.  Coupled  5qll-Yaw  Autopilot 

The  trajectory  sensitivity  functions  A13,/i15,  and 
^18  in  the  linear  case  and  the  correspondent  trajectory 
sensitivity  functions /1175,  /1174,  and/i173  in  the  nonlinear 
case  show  strong  effect  on  the  overshoot  of  the  state  vari¬ 
able  r. 

The  trajectory  sensitivity  functions  <123, /i25,  and 
1.28  in  the  linear  case  and  the  correspondent  trajectory 

sensitivity  functions <1185,  <1184,  and  1.183  in  the  nonlinear 
case  show  strong  effect  on  the  overshoot  of  the  state  vari¬ 
able  p. 

The  trajectory  sensitivity  functions  ^133  ,  135,  and 
138  in  the  linear  case  and  the  correspondent  trajectory 

sensitivity  functions  1165,  ^1164,  and  1163  in  tie  nonlinear 
case  show  strong  affect  on  the  overshoot  of  the  state  vari¬ 
able  ^  . 

The  trajectory  sensitivity  functions  173  ,  175,  and 
178  in  the  Linear  case  and  the  correspondent  trajectory 

sensitivity  functions  1125,  1124,  and  1123  in  the  nonlinear 
case  show  little  effect  on  the  overshoot  of  the  state  vari¬ 
able  So  • 
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The  trajectory  sensitivity  functions  i.83  ,/185,  and 
^88  in  the  linear  case  and  the  correspondent  trajectory 
sensitivity  functions  2.135,  ^1134,  and^l133  in  tae  nonlinear 
case  show  strong  effect  on  the  overshoot  of  the  state  vari¬ 
able  Sd  . 

The  trajectory  sensitivity  functions  /d103  and  .2.105 
in  the  linear  case  and  tha  correspondant  trajectory  sensi¬ 
tivity  functions  /1145  and  ^144  in  the  nonlinear  case  show 
strong  effect  on  the  overshoot  of  the  state  variable  • 
The  trajectory  sensitivity  function  AlOd  in  the  linear  case 
and  the  correspondent  trajectory  sensitivity  fj notion /i143 
in  the  nonlinear  case  show  little  affect  on  the  ever shoot  of 
tha  state  variable  Syc  • 

The  trajectory  sensitivity  functions  /1 113,  2.115, 

and  2.118  in  the  linear  case  and  the  correspondent  trajec¬ 
tory  sensitivity  functions  2- 155,  2.154,  and  2-153  in  the 

nonlinear  case  show  strong  effect  on  the  overshoot  of  the 
state  variable  Sy  • 

The  trajectory  sensitivity  functions  2  123,  2.125, 

and  2128  in  the  linear- case  and  the  correspondent  trajec¬ 
tory  sensitivity  functions  2195,  2.194,  and  2193  in  the 

nonlinear  case  show  strong  effect  on  the  overshoot  of  the 
state  variable  ^  . 

The  trajectory  sensitivity  functions  have  the 
following  range  of  values  : 

Hiniaua  Maximum 

Linear  -  385.25  632.35 

Nonlinear  -  642.9  1772.9 


D.  GEMBBAL  COHCLOSIONS 


The  paraneter 
alaost  all  case. 


strongly 


affect 


the  overshoot  in 
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The  parameter 

case. 

little  affect  the 

overshoot  in 

all 

The  parameter 

etlmost  all  case. 

strongly 

affect 

the 

07 ershoot 

in 

The  parameter 

almost  all  case. 

strongly 

af  feet 

the 

oirershoot 

in 

The  parameter 

strongly 

affect 

the 

overshoot 

in 

almost  all  case. 


E.  BECOHHENO&TIOHS 

It  was  not  possible  to  run  a  CSMP  computer  program  using 
all  stare  variable  (19)  and  parametars  (10)  in  the  nonlinear 
case  due  to  the  worlc  area  available.  The  above  restriction 
occurred  just  when  the  actual  and  nominal  output  were 
required  to  be  printed.  The  computational  time  in  this  case 
may  be  reduced  using  Fortran  subroutines  imbedded  in  the 
CSSP  program  in  order  to  calculate  the  necessary  parameter 
derivatives. 

In  the  Linear  cases  the  computational  time  can  be 
reduced  by  means  of  using  the  method  of  Sensitivity  Points 
that  uses  one  model  for  all  parameters  instead  of  using  one 
model  for  each  parameter. 

The  present  analysis  can  be  repeated  for  the  case  of  the 
Circular*  Airframe  given  in  the  [Hef.  2],  and  comparisons 
between  both  airframes  can  be  performed. 

Future  study  can  be  made  modeling  the  system  in  the 
frequency  domain  where  the  "root  sensitivity"  is  analysed. 

Influence  of  parameter  variations  on  miss  distance  can 
be  analysed  using  an  augumented  system  and  a  scenario  of 
reference. 


*Tha  current,  CSHP  programs  in  th$  appendices  are 
prepared  to  run  the  Circular  case  just  by  inserting  the  data 
of  the  correspondent  airframe. 
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BISSILE  SIZIHQ,  HASS  PBOPBBTIBS  AND  AEBODTMAIIC  DATA 


A.  IMTRODOCTION 

In  this  appendix,  one  considers  a  node!  that  was 
assuned  to  ba  1/6  scale  of  the  actual  elliptical  missile 
configuration  as  given  in  Pig.A.  1. 

Bissile  configuration  was  sized  to  provide  realistc  mass 
properties  needed  for  this  stud;. 

Since  the  main  purpose  of  this  worh  was  to  perform  the 
parameters  sensitivity  analysis  applied  to  a  3 ank-to-Turn 
missile,  no  effort  was  expended  on  a  detailed  design  and 
analysis  of  the  various  parts  envolved,  which  is  given  in 
(Hef.  2].  Some  of  figures  given  in  this  appendix  were  taken 
from  the  (Bef.  2]  fcr  easy  visualization  of  the  system  in 
study. 

B.  GEOHETBIC  AND  HASS  PROPEBTIES  DP  BISSILE  CONPIGDBATION 

One  can  find  in  Table  VI  the  properties  irsed  in  the 
development  of  the  equations  for  applying  the  Parameters 
Sensitivity  Analysis, 

Fig.  A.  2  shows  the  amic  sign  convention,  nomencla¬ 
ture  and  body- fixed  e  the  following  three  assump¬ 

tions  were  made: 

-  Plane  x^  A. 2  is  the  maneuver  plane. 

-  Hissile  ii  in  pitch  (i.  e.  ,  Hy«0  ,  at  fixed 

values  of  ,  and  )  . 

-Hissile  roll  .ate  is  constant. 

Fig.  A.  3  shows  a  block  diagram  of  a  BTT  autopilot. 
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Figara  &.1  Hodal  of  elliptical  coaf igaratioa  (  1/6  scale  ) 


Inertial  acceleraticn  conaands  are  applied  in  pnlar  coordi¬ 
nates  (i.e.,  aagnitode  of  coaaand  applied  to  the 

pitch  autopilot  and  the  direction  (  applied  to  the 

roll  autopilot)  . 

The  yaw  autopilot  is  slewed  to  the  roll  autopilot  to  aini- 
aize  sideslip  angle  by  coordinating  the  aissile  yaw  and  roll 
aotion. 

Achieved  aaneuwer  plane  or  inertial  acceleration  in  rectan¬ 
gular  coordinates  (i.e.,  and  ■^Yw)  is  .deterained  by 
resolving  achieved  body-fixed  accelerations 

(i.e.r**),  and  >)  )  through  aissile  roll  angl8((^  )  (i-e. , 

Euler  angles  e-  and  are  assuaed  to  be  sufficiently 

saall) . 


171 


iiPPENDH  B 
LINE&B  SYSTEM  M30ELS 


A.  IIIRODOCTIOH 

The  following  appendix  addresses  the  informations  of  the 
linear  uncoupled  pitch  channel  autopilot  and  the  roll-yaw 
coupled  channel  autopilot  of  a  elliptical  airframe  taken 
from  CBef.  2]. 

The  linear  tine  domain  analysis  of  the  CBTT  auropilot  was 
assumed  that  the  missile  is  initially  in  the  desired 
maneuver  plane  and  trimmed  at  ten  degrees  angles-of-attack 
(i.e.,  the  eguilibrium  or  trim  angle-of-attack  in  the 

model  of  Fig.A.1  equals  10  degrees  and  the  equilibrium  roll 
rate  Pe  in  the  models  of  Fig. A. 3  agual  zero). Hh an  Pe  =  0, 
Qe  (i.e.,  equilibrium  pitch  rate)  has  been  formed  to  have 
negligible  influence  in  the  lateral  model  compared 
to  o(  and  was  therefore  set  equal  to  zero. 

B.  OHCOOPLED  PITCH  CHAMHBL  AOTOPILOT 

A  general  block  diagram  of  an  uncoupled  pitch  channel 
autopilot  is  shown  in  Fig.  B.1. 

A  normal  acceleration  command  g's)  is  applied  to 

the  pitch  control  law  which  uses  measurements  of  missile 
body  pitch  angular  rate  (q)  and  pitch  normal  acceleration 
(V,)  to  determine  the  required  actuator  command  (&p^)  •  The 
actuator  is  modeled  as  a  first-order  lag  at  3  0Hz  (188.4 
rad/sec) . 

The  dynamic  model  is  linearized  about  a  trim  angle-of- 
attack  as  described  in  appendix  A. 

From  the  block  diagram  of  Fig.  B.1 ,  the  parameters  of 
interest  for  the  present  analysis  are  given  respectively 

**ot  Njp  m^p 
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For  ease  notation  the  correspondence  of  state  variables  and 
constants  are  given  in  Table  VII 

Perforaing  block  diagran  manipulations  see  [Bef.  5],  one 
yields  to  the  following  state  variable  equations: 

Aerodynamic  equations 

11  =  C2  A3  12  ♦  C2  A4  X3  (B.1) 


12  »  XI  -  K  Cl  A1  12  -  K  Cl  A2  X3  (B.2) 

Control  law  equations 

X3  »  -  C3  X3  ♦  C3  Conv  X6  (B.3) 


X4  a  -  Cl  C4  A1  X2  -  Cl  C4  X3  -  C4  X4  '<B.4) 

X5  a  C7  X4  -  C5  X5  -  C6  NZC  (B.5) 

Actuator  equation 

X6  »  C  E  A7  X3  -  6  X6  ♦  C  E  AS  X8  C  E  A6  X11  (B.6) 


C.  COOPLED  BOLL-TAB  CHABXBL  AUTOPILOT 

A  general  block  diagram  of  a  coupled  roll- yaw  channel 
autopilot  is  shown  in  Pig. B.2. 
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The  block 

diagraa  (Fig.B.2) 

sho  ws. 

eight  paramerers  of 

interest  that 

are: 

^  ^  ~  ^ YSy  ' 

A2  =  Qy^  ,  A3  = 

‘5  =  , 

16  .  ,  17  - 

and  A8  =  C  . 

For  ease  of  aotation  the  correspondence  of  state  vari 
ables  and  constants  is  given  in  Table  VIII. 


Eqns.B.7  through  B.17  give  the  state  variable  equations 
of  the  coupled  roll-yaw  channel  autopilot: 

Aerodynanic  aquations 

XI  =  C  Conv(A4  X3  ♦  A3  X8  ♦  A5  X11>  (B.7) 


12  =  C  Conv(A7  X3  ♦  AS  X8  ♦  A6  XII)  (B.8) 

X3  ®  -  XI  -  lalpha/Conv)  X2  (B.9) 

♦  KB  A  A2  X3  ♦  KB  A  Al  XII 

'x4  =  -8  X4  -17.6  X12  ♦  17.6  PHC  (B.  10) 

Control  law  equations  (roll) 

X5  =  -  (0.755/Conv)  X2  -  C  D  A7  X3  (B.  11) 

♦  (0.755  -  8  D)  X4  -  5  X5  -  C  D  A3  X8  -  C  D  A6  XII 
-  17.6  D  X12  +17.6  PHC 
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X6  =  C  E  A7  X3  -  6  X6  ♦  C  E  A8  X8  C  E  A6  XII 


(B.12) 


Z7  =  -  F  KC  (0.755/Conv)  X2  -  (D  E)  F  KC  A7  X3  (B.  13) 

♦  P  KC(0.755  -  8  D)  X4  ♦  KC(15  -  5  F)  X5 

♦  KC(6P  -  15)  X6  -  15  X7  -  (D  ♦  E»  P  KC  C  A8  X3 
”  (D  ♦  E)  P  KC  C  A6  XII  -  P  KC  D  17.6  X12 

♦  P  KC  D  17.6  PHC 

X8  =  188.4  Conv  X7  -  188.4  X8  (B.14) 

Control  law  equations  (yaw) 

X9  =  (K1  A  A2/T1)  X3  -  X9/T1  ♦  (K1  A  a1/T1)  XII  (B.15) 


XI 0  *  (K2/Conv)  XI  -  (K2  H  ALPHAS) /Conv  X2  (B.  16) 

♦  (K2/tO)((Kl  A  A2)/T1  ♦  B  A4  -  3  ALPHAB  C  A7  X3 

♦  (K2/10)  (B  A3  -  e  ALPHAB  C  A8  X3  ♦  K2(1  -  (VIO  T1)  X9 

♦  (K2/10)  ((K1  A  A1)/T1)  ♦  B  A5  -  H  ALPHAB  C  A6)  XII 

Actuator  equations 

XII  =  188.4  conv  X10  -  188.4  XII  (B.  17) 


XI  2  =  X2/Conv  (B.  18) 
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D.  AEBODTIiHIC  D&ri-LINEiB  APPEOXIE ATIOH 


&  linear  approach  was  used  in  the  design  and  stability 
analysis  of  the  aultopllots  of  tha  pitch,  yaw,  and  roll 
channels,  both  uncoupled  and  coupled.  According,  a  linear 
approximation  of  the  aerodynamic  darivatives  at  m=3.95  was 
provided,  about  which  the  system  could  be  perturbed.  These 
linearized  aerodynamic  derivatives  are  given  in  Q  Eef .  2]  and 
are  presented  here  in  Table  ix. 


T&BLB  TII 

Correspondence  of  Syebols  (Onconpled  Pitch  hatopilot) 


q 

«Sp 

X 

Y 

<5  Pc 

«.p 

qs/d 

57,3qsd/I 

188.4 

150.0 

6.0 

0.53544 

0.48 

0.38375 

3.07 

0.48 

57.3 


XI 

X2 

X3 

X4 

x5 

X6 

k^ 

&2 

&3 

kU 

Cl 

C2 

C3 

C4 

C5 

C6 

C7 

C8 

C9 

K 

ConT 
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TABLE  Till 

Correspondence  of  Symbols  (Complad  Boll*Taw  Aatopilot) 


C 


cV 

Hr 

%y 


3.  1416 
57.3 
qs/d 

qs^L^/l22 
qsdl/lKx 
0.05033 
0.078 
0.25 
oc  e 


Conv 


0.25 

0.84 

6.08 


4.17 

0,48 

2,0 


ALPHA 


oCg  /57.3  .  ALPHAS 

0.48  .  KB 


APPEMDII  C 

NOHLIHEAR  STSTEH  MODEL 


A.  IITBODOCTIOH 

A  general  blocic  diagraa  o£  the  nonlinear  noiel  is  shown 
in  Fig. C.1  The  three  dimensional  nonlinear  aerod/namic  model 
presented  here  is  for  the  same  conditions  used  in  the  linear 
model  shown  in  appendix  B.  The  same  flight  condition  ased 
for  the  linear  case  is  the  same  for  the  nonlinear  model  (i. 
e.  ,  60  kft  altitude,  Mach  3.95).  The  control  laws  are  the 
same  used  for  the  linear  models  except  for  a  minor  modifica* 
tion  to  the  coordinating  branch  dependence  on  angle-of- 
attack  and  also  the  inclusion  of  anti-gravity  bias  as  stated 
in  [Ref.  2]. 


B.  COITBOL  LABS 

The  control  laws  used  for  the  following  nonlinear  3-D 
studies  were  the  same  as  those  used  for  the  linear  models, 
except  for  the  gain  X  shown  in  the  bold  line  of  coordina¬ 
tion  branch  in  Pig.  C.1.  The  new  gain  ^  is  held  constant 
at  one  degree  magnitude  for  angle-of-attack  less  than  one 
degree  positive.  For  angle-of-attack  greater  than  one 
degree  positive,  the  gain  oc  is  equal  to  the  angle-of- 
attack.  This  keeps  coordination  for  very  small  angle-of- 
attack. 

Gravity  effects  were  not  included  in  the  linear  models 
in  appendix  B,  because  it  was  assumed  to  have  a  negligible 
influence  on  autopilot  stability  and  response  for  pertuba- 
tions  about  a  missile  trim  condition.  However,  gravity 
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Figure  C. 1  lonllne&c  Hodel 
185 


effects  were  included  in  the  following  nonlinear  model  where 
the  missile  body-fixed  yaw  axis  will  be  subjected  to  the 
full  force  of  gravity  and  may,  therefore,  have  a  significant 
influence  on  sideslip. 

The  gravity  effects  are  compensated  for  by  pitch  and 
roll  notion  of  th^  missile  which  have  less  influence  on 
sideslip  than  yaw  motion.  In  inertial  rectangalar  coordi¬ 
nates  one  has: 

^acceleration  command  in  inertial  direction 


=  ^7  -  Cos(e) 

‘C 


(C.1) 


where 

=  guidance  connand(  g*s  ) 

-Cos®  =  anti-gravity  bias  command  (  g's  ) 

=  acceleration  command  in  the  inertial  direction 
guidance  command (  g*s  ) 

There  is  no  gravity  effect  in  I  direction.  In  order  to 
assure  an  anti-gravity  bias  of  just  one  g,  oa  e  needs  to 
modify  the  anti-gravity  command  as  follows: 

elliptical  airframe  inertial  acceleration  command 

~  Nzc-0. 91  3cos  (()) 

C.  PEBFOBAHCB  COBHiHOS 

The  commands  2  gees  is  first  applied  in  the  Qo  or 
upward  direction  at  2  seconds.  Since  both  the  missile 
roll  angle  and  roll  angle  command  are  at  zero  degrees,  there 
is  no  roll  motion  and  the  missile  turns  upward  as  a 
shid- to-turn  controlled  missile.  &t  5  seconds,  a  second 
2  g's  inertial  guidance  command  is  applied  in  the  downward 
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-r 


or  1800  direction.  The  missile  is  commanded  to  coll  through 
1800  while  mowing  in  a  coordinated  manner  in  yaw  and  roll  to 
minimize  sideslip  angle  and  prevent  or  minimize  negative 
angle -of -attack. 

The  response  of  achieved  maneuver  and  cross-plane  accel¬ 
erations  during  the  first  2  seconds  are  due  to  initial 
conditions,  gravity  and  anti-gravity  bias  effects.  The 
initial  conditions  were  added  to  minimize  the  transients 
which  result  when  gravity  bias  commands  the  autopilot  for 
constant  alritude  missile  flight.  The  initial  conditions 
were : 

=  angle- of -attack  =  2.41  deg 
Sp  3  pitch  tail  angle  =  0.658  deg 
©  =  pitch  Euler  angle  »  3.65  deg 

output  of  pitch  acceleration  feedback  lag  =  -1.0  g's 
pitch  control  law  acceleration  error 
lag  prior  to  dc  gain  =  -  0.3105 
=  pitch  actuator  command  =  0.S58  deg 
The  achieved  maneuver  plane  acceleration  (Bgn.C. 2)  is 
calculated  from  the  body-fixed  acceleration 
the  roll  angle  {)  as  follows: 


’l.a'  \ 


6 


and 


During  the  first  command,  achieved  body-fixed  yaw  accel¬ 
eration  (  and  missile  roll  angle  are  equal  to  zero, 
because  the  roll  channel  is  not  commanded.  Therefore, 
achieved  maneuver  plane  acceleration  is  equal  to  the  body- 
fixed  acceleration  'V\  z^.  During  the  second  command  the 
missile  roll  angle  has  the  same  roll  angle  response  of  the 
roll  autopilot. 
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0.  &EB0DTIiHIC  D&TA  -  HOMLIMB&B  BEPBES EMTATIOH 


The  aerodynaoic  data  are  given  in  [Bef.  2].  The  entire 
study  was  conducted  using  t!=»3.95  and  the  aerodynamic  coeffi¬ 
cients  are  based  on  a  body- fixed  axis  system  of  Pig.B.2. 

For  reference,  the  normal  force  and  pitching  moment 
plots  given  in  appendix  B  of  [Bef.  2]  at  8=3.  95  were  repro¬ 
duced  in  Table  X  and  XI,  which  give  just  the  values  for 
points  of  interest.  These  tables  present,  too,  the  corre¬ 
spondent  derivatives  that  were  found  by  simulation  using  a 
fortran  program  shown  in  appendix  P.  The  aerodyiamic  deriv¬ 
ative  of  Cjj  ,  ,  and  Cy  with  respect  to  sideslip 

auigle  (h  ,  yaw  control  6y  ,  and  roll  control  aure 

presented  in  Table  XII,  XIII  and  XIV  which  give  just  values 

for  points  of  interest  and  includes  the  correspondent  deriv¬ 
atives  that  ware  found  by  simulation  using  a  fortran  program 
given  in  appendix  P. 

The  block  diagram  of  Pig.C.1  shows  10  parameter  of 
interest  that  are  ; 

A1  =  C^(«c,fp)  ,  A2  =  C^(x,*p)  ,  A3  =  , 

A4  =  C  (oc)  ,  A5  =  C  («.)  ,  A6  *  Cy  (  <x. )  , 


=  Cc^^(o^)  ,  AS  =  Co,  (X)  ,  A9  =  (o<  )  , 


and  AA  =  Cy.  («.)  . 


Osing  the  same  procedure  as  given  in  appendix  B  one  can 


have  the  Tables  X7  and  XVI  that  give  the  correspondence  of 
symbols  for  the  nonlinear  model.  Bgns.C.3  through  C.21  give 
the  state  variable  equations  of  the  nonlinear  model. 
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XI  «  -  Ctt  XI  -  C02  C4  CN(X6,X4) 


(C.3) 


X2  *  C7  XI  -  C5  X2  ♦  C6  C19  Cos(Xl7)  -  C6  NZC  (C.4) 


X3  =  -  C7  C8  XI  ♦  (C5  C8  -  C9)  X2  ♦  (C9/Conv)  X5  (C.5) 


♦  C6  C8  NZC 

X4  *  -  C3  X4  ♦  C3  Conv  X3  {C.6) 

« 

X5  =  (X17  X18)/Conv  ♦  Cl  Conv  CM(X6,X4)  (C.7) 

X6  =  X5  -  KB  CH(X6,X4)  -  (X16  X18|/ConT  (C.8) 

X7  =  X5  Cos(X19) /Conv  -  X17  Sin  (XI  9 ) /Conv  (C.9) 


X8  *  -  C10  X8  ♦  (C12  -  C11  C13)  X10 


(C.  10) 


♦  C11  C14  PRC  -  C11  C14  X19  -  CO  Cl  l(Cj,^  (X6)  XI  5 

♦  Co,  (x6)  XI 3)-  C12  X18/Conv  ^ 


X9  »  -  X9/T1  ♦  (KI/TI)  C02(C^  (X6)  X16 

^  0 


(C.  11) 


♦  C^  (X6)  X15) 

Tffy 
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X10 


C13  no  ♦  C14  PHC 


C14  X19 


(C. 12) 


XII  =  -  C15  XII  ♦  C16  C0(Cj^(x6)  X16  (C.13) 

♦  C  (X6)  X15  ♦  C  (16)  X13) 

X12  -  -  Cl  X12  ♦  K  (C17  -  C10/C18)  X18  (C.  14) 

♦  (K/C18)(-12  -  C11  C13)  X10  ♦  K  C11  (C14/C18|  PHC 

-  K  C11  C14/C18  X19  -  (K/C18)  CO  C  (X6)(C11 

♦  C16)  X16  -  (K/C18)  CO  Cj^(X6)(C11  ♦  C16)  X15 

-  (K/C18)  CO  Cj,^(X6)  (Cl  1  ♦  C16»  X13 

-  (K  C12/C3nv)  X18  ♦  (C15/C18  -  C17)  111 

*X13  =  -  C3X13  ♦  C3  Conv  X12  (C.  15) 


X14  =  (K2/10)  (-  X9/T1  ♦  K1  T1  C02  (C  (X6)  X16  (C.16) 

♦  Cy^^(X6)  X15)  -  fX15  XI  8) /Conv2  (X5  -  KB  C02  CN{X6,X4) 

-  (X16  X18)/Cor  i18  -  (KYP/Con\r)  (X5  -  KB  C02  Cm(X6,X4) 

-  (X16  X18)/Conv  X18  -  (KlP/Conv)  CO  XX6  (C^^(X6)  X16 

-  Cjj^(X6)  115  ♦  C,j^(X6)  X13)  ♦  (K2/Conv)  X 17 

♦  K2  X9  -  (K2  KYP/Conv2)  XX6  X18 

# 

X15  =  -  C3  X15  ♦  C3  conv  X14  (C.  17) 


X16  »  KB  C02(C  (X6)  X16  ♦  C^  (X6)  X15 

♦  (X16  X18)  /Conv  -  X17 


(C. 18) 
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X17 


(X5  Xl8)/Ccnv  ♦  C01  ConT(3^ 


(X6) 


X16 


♦  (X6)  X13) 

X18  =  CO  coav  (Cj^^  (X6)  X16  ♦ 

♦  C  (X6)  X13) 

XI 9  =  X18/C3nv 


X15 


(C.  19) 


(C.  20) 


(C.  21) 
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microcopy  resolution  test  chart 
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